MASSIMILIANO LOCATELLI 
Mid infrared digital 


holography and 
terahertz imaging 


DIOXTREKTOIRTASIKO 
VERSIT, SLI STUDI DI FIRENZE 


= 


PARSESMIEO anzi 
FIRENZ NIVERSITY P 


iy N 
ì\ 
c A 
Z 
| 
IS 


| 


ANY, 
gS 
DA 
FIRENZE 


PRESS 


PREMIO TESI DI DOTTORATO 


— 42- 


PREMIO TESI DI DOTTORATO 


Commissione giudicatrice, anno 2013 


Luigi Lotti, presidente della Commissione 


Tito Arecchi, Area Scientifica 

Franco Cambi, Area Umanistica 

Paolo Felli, Area Tecnologica 

Michele Arcangelo Feo, Area Umanistica 
Roberto Genesio, Area Tecnologica 

Luigi Lotti, Area Scienze Sociali 

Mario Pio Marzocchi, Area Scientifica 
Adolfo Pazzagli, Area Biomedica 

Mario Giuseppe Rossi, Area Umanistica 
Salvatore Ruggieri, Area Biomedica 
Saulo Sirigatti, Area Biomedica 

Piero Tani, Area Scienze Sociali 

Fiorenzo Cesare Ugolini, Area Tecnologica 
Vincenzo Varano, Area Scienze Sociali 
Graziella Vescovini, Area Umanistica 


Massimiliano Locatelli 


Mid infrared digital holography 
and terahertz imaging 


Firenze University Press 
2014 


Mid infrared digital holography and terahertz imaging 
/ Massimiliano Locatelli. - Firenze : Firenze University 
Press, 2014. 

(Premio Tesi di Dottorato; 42) 


http://digital.casalini.it/9788866556893 


ISBN 978-88-6655-688-6 (print) 
ISBN 978-88-6655-689-3 (online) 


Peer Review Process 

All publications are submitted to an external refereeing process under the responsibility of the FUP Editorial Board 
and the Scientific Committees of the individual series. The works published in the FUP catalogue are evaluated and 
approved by the Editorial Board of the publishing house. For a more detailed description of the refereeing process we 
refer to the official documents published in the online catalogue of the FUP (www.fupress.com). 


Firenze University Press Editorial Board 


G. Nigro (Co-ordinator), M.T. Bartoli, M. Boddi, R. Casalbuoni, C. Ciappei, R. Del Punta, A. Dolfi, V. Fargion, S. 
Ferrone, M. Garzaniti, P. Guarnieri, A. Mariani, M. Marini, A. Novelli, M. Verga, A. Zorzi. 


This work is licensed under a Creative Commons Attribution 4.0 International License (CC BY 4.0: 
http://creativecommons.org/licenses/by/4.0/) 


CC 2014 Firenze University Press 
Università degli Studi di Firenze 
Firenze University Press 

Borgo Albizi, 28, 50122 Firenze, Italy 
www.fupress.com 

Printed in Italy 


Table of contents 


Introduction 


Chapter 1 
Analog and digital holography 


1. Analog holography 

1.1 Wavefront recording 

1.2. Wavefront reconstruction 

1.3 On-axis and off-axis holography 
2. Digital holography 

2.1 Wavefront recording 

2.2 Wavefront reconstruction 


2.3 Advantages of digital holography 


Chapter 2 
Infrared digital holography 


1. Infrared radiation, sources and detectors 
1.1 Infrared radiation 
1.2 Sources 
1.3 Detectors 

2. Advantages and disadvantages of IR DH 


Chapter 3 
Mid IR digital holography experiments 


1. Experimental apparatus 
2. Experimental setups and results 


2.1. First results 


2.2 Super resolution by means of synthetic aperture 


2.3 Human size holograms 


2.4 Visible analog reconstruction of IR digital holograms 


2.5 Smoke and flames hidden object holograms 


Massimiliano Locatelli, Mid infrared digital holography and terahertz imaging, ISBN 978-88-6655-688-6 (print) 


ISBN 978-88-6655-689-3 (online), CC BY 4.0, 2014 Firenze University Press 


Mid infrared digital holography and terahertz imaging 


Chapter 4 
Terahertz imaging experiments 


1. Experimental apparatus 

2 Experimental setups and results 
2.1 THz imaging with single point detector 
2.2 Toward Terahertz digital holography 


Chapter 5 
Conclusions 


List of Publications 


Bibliography 


65 
65 
67 
67 
76 


79 


81 


83 


Introduction 


This research project was devoted to the development of imaging techniques in 
the Mid Infrared and Terahertz regions. As we will see, these portions of the elec- 
tromagnetic spectrum offer interesting opportunities for imaging purposes under 
various points of view. In particular we focalized our attention toward the Mid Infra- 
red radiation generated by a high power CO; laser working at 10.6 um and toward 
the 10 times longer Far Infrared radiation generated by a Quantum Cascade Laser 
working at 107 um. As we will see, each of these two long wavelength radiation 
sources has peculiar characteristics that can be fruitfully exploited for imaging. In 
this work, in particular, we decided to explore the possibilities offered by these 
sources in a specific imaging technique, Digital Holography. This is a quite recent 
technique (its birth can be dated back to 1994 by Schnars and Jüptner) and Infrared 
radiation digital holography is an even more recent research field. In particular, the 
first digital holography experiments at 10.6 um were performed at the National In- 
stitute of Optics (INO), in Florence, about in 2003 but it is only very recently that 
this technique has started receiving new attention because of its possible application 
to real time non-destructive testing of large size samples and because of its possible 
involvement in many important applicative fields like homeland security and real 3d 
television. During this research project digital holography in the Mid Infrared has 
been deeply investigated and innovative configurations and applications have been 
developed. Digital holography in the Terahertz region is an almost unexplored re- 
search field, the first results being achieved very recently and under limitative condi- 
tions but, as every Terahertz imaging system, it is rich in potential applications. Dur- 
ing this research work we implemented a simple standard imaging technique in this 
spectral region and we used it to investigate the performances of a new kind of Te- 
rahertz radiation detectors: nanowire detectors. Finally, we traced the research lines 
we intend to develop in the field of Terahertz digital holography using QCLs radia- 
tion and microbolometric focal plane arrays. 

The experiments presented in what follows were obtained in the framework of 
collaboration between the European Laboratory for Non Linear Spectroscopy (Lens), 
the CNR National Institute of Optics (INO-CNR), and the CNR Nanoscience Insti- 
tute of Pisa (CNR-NANO), who provided the Terahertz source and the innovative 
detectors tested in this work. 

In the first chapter we briefly review the basic principles of analog holography 
and describe, in more detail, its natural development, digital holography. In particu- 
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lar, in the last part of the chapter, we discuss the peculiar features that characterize 
digital holography with respect to the other imaging techniques. 

In the second chapter we introduce the basic notions about the two Infrared re- 
gions employed in our investigations and give a general description of some of the 
most interesting sources and detectors working in each of these two spectral ranges. 
We then focus on the possibility and opportunity to use Infrared radiation to per- 
form digital holography and discuss the advantages, and disadvantages, that this 
spectral region offers in this field. 

In the third chapter we describe and discuss the experimental apparatus and re- 
sults obtained, during this research project, in the field of Infrared digital holography 
at 10.6 um. In particular, in the paragraph about the experimental results, five differ- 
ent experiments are illustrated in detail. The last of these experiments deals with a 
novel and promising capacity of Infrared digital holography, discovered during our 
tests: the possibility to see through smoke and flames. This result, as can be easily 
understood, could have important applications in the homeland security field (to 
help Firefighters in fire scenarios) as well as in the industrial field (to perform non- 
destructive tests during high temperature and flames involving procedures) and has 
stimulated the attention of the scientific community toward Infrared digital holog- 
raphy. 

In the fourth chapter we describe and discuss the experimental apparatus and re- 
sults obtained, during this research project, in the field of Terahertz imaging; in par- 
ticular we illustrate the results obtained with a simple standard imaging system and 
examine the possibility, and the opportunity, to translate the experimental configu- 
rations developed in the field of Mid Infrared digital holography to the Terahertz 
range. 

In the fifth chapter we summarize the results obtained during this research pro- 
ject, we discuss the limitations of the investigated techniques and discuss their possi- 
ble future developments and their potential applications. 


Chapter 1 
Analog and digital holography 


Holography was ideated in 1948 by the naturalized British Hungarian physicist 
Dénes Gabor, during his studies on the electronic microscope [1]. Gabor devised a 
technique to record and, subsequently, reconstruct both the amplitude and the phase 
information of a light wavefront diffracted by an object irradiated with coherent ra- 
diation; Gabor understood that this remarkable result could be obtained simply ex- 
ploiting the interference pattern resulting from the superposition of the radiation 
coming from the object (object beam) with a portion of the radiation used for irradi- 
ating it (reference beam). Because of its ability to record the image of the object as a 
whole, the technique, initially called by Gabor simply “Wavefront Reconstruction’, 
acquired the name of ‘Holography’, from the Greek words ‘holos’, meaning whole, 
and ‘graphein’, meaning writing. The technique devised by Gabor was greeted ini- 
tially with little interest by the scientific community, but starting by the '60s, that is 
when high spatial and temporal coherence radiation became easily accessible thanks 
to the invention of the laser, it underwent an important theoretical and experimental 
development so that Gabor, in 1971, was awarded the Nobel Prize in Physics for his 
invention. 

In analog, or classical, holography, the interference pattern obtained by the su- 
perposition of the two beams is usually recorded on an analog image recording de- 
vice (usually some appropriate kind of photographic film) that, once developed, con- 
stitutes the so-called ‘hologram’; the object is then removed from its position and the 
transparency is repositioned at the point where it had been impressed; if the photo- 
graphic plate is again irradiated by the reference beam, one observes the reconstruc- 
tion, for diffraction, of a wavefront identical to the one coming from the object be- 
fore it was removed; this means that the observer is invested by a light wavefront in- 
distinguishable, with passive means, from the one generated by the object itself, and 
then, within a certain viewing angle, equipped with the same perspective and depth 
of focus. 

Technological developments in computing and in image recording devices, have 
soon suggested the idea of transferring the analog recording and reconstruction pro- 
cesses of classical holography into the digital domain, giving rise to new techniques 
such as Computer Generated Holography, in which artificial interferograms of non- 
real objects, numerically created, are used to realize real holograms to be analogically 
reconstructed, and Digital Holography, in which the recording process is carried out 
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by means of digital recording devices and the reconstruction process is performed 
numerically. The first step towards the realization of the latter technique was made 
in 1967 by Goodman and Lawrence [2] who managed to numerically reconstruct, in 
about 5 minutes, a digital hologram obtained from a real hologram recorded on a 
photographic film. The decisive step towards the complete digitalization of the holo- 
graphic process was carried out in 1994 by Schnars and Jiiptner [3] which realized, 
for the first time, the numerical reconstruction of a CCD recorded hologram. 
Thanks to technical progresses of modern CCD and CMOS devices and thanks to 
the higher computing power of today's computers, digital holography is experienc- 
ing massive growth and, because of its methodological simplicity, is replacing, in 
most applications, the corresponding classical technique; this process is further 
speeded up by the fact that the major manufacturers of the most popular photo- 
graphic films used in analog holography have closed down their production lines 
due to the general tendency to prefer modern digital technologies. However, to date, 
this technique still suffers from some limitations due to the lower resolution of elec- 
tronic recording devices compared to analog ones, limitations intended to diminish 
with the advent of bigger and bigger CCDs with increasing resolution. The long time 
required for the numerical reconstruction that, in the early days of this technique, 
constituted a significant limitation for the development of digital holography, have 
now been greatly reduced and modern calculators allow, in many cases, a recon- 
struction substantially simultaneous with the hologram recording procedure. Real 
time wavefront reconstruction, made it possible by digital holography, represents 
today a very important and attractive tool in numerous applications. 


1. Analog holography 


As already mentioned in the previous section, the goal of holography is the re- 
cording and subsequent reconstruction of the amplitude and phase information of 
the wavefront coming from an object irradiated with coherent radiation; holography 
is, therefore, a two-step process: a first step to record the interference pattern, the 
hologram, and a second step to reconstruct the object wavefront or, as is usually 
said, to reconstruct the hologram. 


1.1 Wavefront recording 


Image recording devices provides a response proportional to the average energy 
impinging on their surface; their response is therefore proportional to the irradiance, 
or average intensity, I, of the electromagnetic field; this physical quantity, in fact, 
represents the average energy transported from the field, in the direction of propaga- 
tion of the field itself, divided for the time and surface units; the average is per- 
formed on a time Tm much greater than the period of the electromagnetic wave and 
therefore the irradiance measurement results directly proportional to the square 
amplitude of the field but any information on its phase distribution is lost; for this 
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reason, if the wavefront phase information is of interest too, it is necessary, some- 
how, to bring back the phase measurement to an irradiance measurement. The idea 
at the basis of holography is precisely to exploit the phenomenon of interference be- 
tween electromagnetic waves to incorporate the phase information in the intensity 
information. In holography, in particular, the phase information of the object wave- 
front is ‘written’ in the interference pattern created by the object beam and an ap- 
propriate reference beam. 

Depending on the nature of the object under investigation two different holo- 
graphic configurations, with respect to the origin of the object beam, are possible: 
transmission holographic configuration and reflection holographic configuration. 

For transparent samples the coherent radiation passes through the sample so 
that the ballistic radiation transmitted through the object itself, modified by its 
transparency characteristics, constitutes the object beam (figure 1.1). 


Figure 1.1: General transmission holographic configuration. 


In this case, both the local absorption coefficient of the material which consti- 
tutes the sample and the optical path that the beam travels within the sample itself, 
affect the wavefront emerging from the sample; in the first case, the beam undergoes 
a change in its amplitude while, in the second case, it undergoes a phase change de- 
pending on the local refractive index of the material and on the local thickness of the 
sample. 

For opaque samples the coherent radiation is directed toward the sample so that 
the reflected, or scattered, radiation in the direction of the recording device consti- 
tutes the object beam (figure 1.2). 
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Figure 1.2: General reflection holographic configuration. 


In the ideal case of a totally reflecting object (that is with surface roughness 
much shorter than the wavelength in use) the recording device receives only the ra- 
diation directly reflected from the object according to the laws of geometrical reflec- 
tion; in this case it is possible to record in the hologram only the information of that 
portion of the object wavefront directly reflected toward the recording surface. If in- 
stead the surface of the object presents a roughness comparable to the wavelength in 
use, one can think the object as consisting of a huge number of small reflective ele- 
ments with random orientation each of which irradiate in a different direction and 
one speaks, in that case, of scattered radiation; the various waves emitted by these 
microscopic sources interfere with each other and what is observed on the recording 
surface is a dotted image composed of a set of bright and non-luminous randomly 
distributed small spots, the so called speckles [4]; in principle, any portion of the re- 
cording device can thus collect the information coming from each portion of the ob- 
ject, regardless of its size; in this configuration, called speckle holography, therefore, it 
is possible to use any portion, however small, of the hologram, to reconstruct the en- 
tire wavefront, even though, as we will see, the smaller the portion of the hologram, 
the lower will be the final reconstruction resolution. We have to consider also that, 
in holography, speckles contain the information on the desired object wave but re- 
duce, at the same time, the final reconstruction resolution; in particular, the effect of 
this noise is linked to the size of the speckles observed in the reconstructed image 
plane, namely the so-called subjective speckles, so that the larger the speckles the 
more disturbing they are. It can be shown that a larger subjective speckle size is ob- 
tained by reducing the hologram aperture [4]. Real objects never act as ideal reflec- 
tors and thus, in practice, one can always think of being in a case, more or less 
marked, of scattering object; what changes is just the percentage of scattered radia- 
tion that ends up in the recording device. 
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The object beam and the reference beam usually travel different paths but, in or- 
der to obtain interference fringes, they must be coherent with each other and there- 
fore, usually, are extracted from a single laser beam by means of suitable beam split- 
ters. For the same reason, the two paths must be equal within a value lower than the 
laser coherence length. Furthermore, as in every interferometry experiment, in order 
to get good interference fringe visibility, the two beams should have comparable in- 
tensity and the differences between the two optical paths should be kept stable, with- 
in a fraction of the wavelength, during the time required to perform the hologram 
recording. Finally, as we will see later, in order to obtain an undistorted reconstruct- 
ed object wavefront during the reconstruction step, the reference beam amplitude 
has to be uniform across the recording device surface; strictly speaking every homo- 
geneous reference wave could therefore be used; usually, however, only plane waves 
or large curvature spherical waves are employed. 


If we consider monochromatic waves of frequency w and if we denote with 
E,(x,y) = E(x, ye'o%”) and E,(x,y) = E,e!9r& the complex amplitudes of 
the object and reference beam, respectively, on the generic plane (X,Y), the total ir- 
radiance of their interference field on the recording plane, or hologram plane, 
(Xn Yp), is 


I(Xy, Yy) X [Eo (Xn Yy) + Er Om Yy)? = 
T [Eo (Xp; Ya)? + |E, (Vx) I? + ES (Xu, Yu) Er Xu Yy) + Eo Xu, Yy) Er Cy Yu) = 
= ES (Xn, Vu) + E? + 2E, (xy, Yn )E, cosld, (xu, Yu) — $o (Xn Vu)] 


The first two terms of this last expression depend only on the amplitude of the two 
fields, the third term, instead, depends directly on their relative phase; the interfer- 
ence pattern recorded on the hologram results thus to be a function of both the am- 
plitude and the phase of the object field. The recording medium, whatever it is (sil- 
ver halide photographic slab, photorefractive material, photopolymer, dichromatic 
jelly, etc.), must be used within the linearity region of its transmittance vs exposure 
curve and provides, in this case, a linear mapping of the resulting field intensity, in 
the sense that, once developed, locally changes its physical characteristics (absorp- 
tion coefficient, refractive index, thickness, ...) so as to present a transparency 
H (xy, Yp) linearly proportional to the incident irradiance 


H (xy, Yu) = ho + Bt (xy, Yy) = 
= ho + BILE (xy, Vu)? + |E- (x YH) I? + ES (Xu, Yu) E (xy Yu) 
+ Eo Xn Yu) EF OH Yu)] 
Where 


ho is the transparency of the unexposed plate (not present in digital holography), 
T is the exposure time, 


13 


Mid infrared digital holography and terahertz imaging 


B is a proportionality constant characteristic of the recording medium (negative for 
photographic emulsions). 


1.2. Wavefront reconstruction 


The second step in a holography experiment is the reconstruction of the object 
wavefront ‘recorded’ in the hologram. Once the holographic plate is chemically de- 
veloped it is rearranged in its original position while the object is removed; at this 
point the plate can be considered as a diaphragm equipped with the transparency, or 
transmittance (in general complex), calculated above. The holographic reconstruc- 
tion is then obtained by sending an appropriate beam, the so-called reconstruction 
beam, in order to create, for diffraction through the hologram, a wavefront identical 
to the one previously coming from the object. 

If we send on the transparency a generic reconstruction beam R(x, y), the ex- 
pression of the wave at the exit of the transparency is 


E (xy Vu) = R(xy Yy) A (XH, Yy) = 
= Rp Vy ho + Rw VH)BTE (XH, Vw)? + ROW VWBTE, i VWI? + 
+ R(Xy Yy) PTE (Xu, Yy)Er (XH Vu) + Ry Vo )BTEo Xy Va Er XH Yy) 


If we use, as a reconstruction beam, a beam exactly identical to the one used to 
record the hologram, namely R(x, y) = €,(x,y), the wavefront at the exit of the hol- 
ogram becomes 


EX Yn) = 
= E (xy Yy)ho + Er (Xy, YH)BTE0 (Xn Yy)? + E Cw Va PTE, (x Yu? + 
+E, (xy Yu)BTES Xn Yn)Er Cy Yu) + Er yp Yn)BTE (Xp, VWEr (x Ya) = 
= {[ho + BLES (xu, yu) + E21]}E Gy Yn) + 
+PTEF (x, Ya)E XH Vu) + BLES Eo (Xy, Yy) 


This equality in the plane of the hologram is valid, thanks to the principle of exist- 
ence and uniqueness of differential equation with fixed contour condition (on the 
hologram plane output), in every plane (X,Y) of the observation half-space and we 
can therefore write, more generally 


E(x,y) = {ho + BTIEZ(x,y) + EEFDE, (x,y) + BLE? (x, YEG (x, y) + BTEFE, (x,y) 


The first term of this expression is the so-called undiffracted (or zero-order or DC) 
term and coincides with the reconstruction (and reference) beam E, (x, y) except for 
a multiplicative factor in general dependent on the spatial variables. The interpreta- 
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tion of the second term is a bit more complex but it can be shown [4;5], by means of 
the so-called imaging equations, that it converges to form a real pseudoscopic and 
distorted object image, placed, with respect to the hologram plane, in the half-space 
opposite to the one where the object was situated originally; this image is distorted 
because of the complex factor, depending on the space variables, €? (x, y). The third 
term coincides with the object beam €,(x, y) except for a real multiplication factor 
which influences exclusively, locally, its brightness; this term is just what we needed, 
namely a wavefront identical, except for a uniform multiplication factor, to the one 
coming from the object during the recording process; an observer positioned in the 
half-space opposite to the object original position with respect to the hologram 
plane, can thus see a virtual image, originating beyond the hologram, substantially 
identical to what he perceived in the presence of the ‘real object’, completely indis- 
tinguishable, apart from a possible brightness uniform attenuation, from the original 
one. 

If an undistorted real image of the object is desired, it is possible to illuminate 
the hologram with a reconstruction beam equal to the complex conjugate of the ref- 
erence beam, R(x, y) = E;(x,y), in which case we have 


E(xy, Yu) = [ho + Bi(E? + EP )]E} (xy Yu) + BrEzes (Xn Yn) 
+ BTE}? (xu, Yn) Eo (Xn, Yy) 


Here the undistorted term is just the one relative to E (xy, Yg) which, as we have 
said, produces the real image of the object under examination. 


1.3 On-axis and off-axis holography 


In the configuration originally designed by Gabor, subsequently defined on-axis 
holography, the object and the reference beams were both coaxial with the axis of the 
recording plate; in particular, in the experiment by Gabor, the object was a highly 
transparent object with weak transmittance variations with respect to the medium 
transmittance; in these conditions both the reference and the object beam could be 
directly obtained from the same laser beam which, passing through the sample, was 
split into two parts: a major portion of the beam which continued undisturbed and 
worked as a reference beam and a minor portion of it which was diffracted by the 
transmittance changes of the sample and which constituted the object beam. This 
kind of configuration, however, presents a problem, the so-called twin image prob- 
lem: the three diffraction orders forming the reconstructed wavefront (the real im- 
age, the virtual image and the undiffracted order) are spatially superimposed in the 
observation direction, troubling each other. To overcome this drawback, in 1962 
Leith and Upatnieks devised what was defined off-axis configuration [6]; in this con- 
figuration the reference beam and the object beam are not coaxial with the axis of 
the recording plane but are inclined to each other by an angle @; with the increasing 
of this angle the interference fringes frequency increases rapidly until fringes become 
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indistinguishable to the naked eye, and, as we shall see, in the hologram reconstruc- 
tion, the real image, the virtual image and the zero order term spatially separate pro- 
gressively in proportion to this inclination angle. Assuming that the object is posi- 
tioned in front of the recording plate so that the object beam is mainly directed along 
the hologram axis (Z axis) and assuming that the reference beam is a plane wave in- 
clined with respect to the hologram plane by an angle @ along the Y axis, (k = 
ksin0$9 +kcos@ 2), we can write, on the hologram plane (z = 0), 


En y) = Bey 


sin(@) 


with a = A 


The intensity pattern resulting from the interference of the two beams on the 
hologram plane is therefore 


1x4 Yn) = 
= E, (Xp, Vu? + [Epe2" |” + E; (4, Yn)EreT +E, (xy, yy )Ee eR = 


= E? (xu, Yy) + E? + 2E, (Xy, Yy )Er cos[2ray — po (Xy, Vu) 


This expression shows that the amplitude and the phase of the wavefront coming 
from the object are recorded, respectively, as amplitude and phase modulations of a 
carrier frequency a@ and, as we shall see, if the carrier frequency is sufficiently high, 
the amplitude and phase distributions can be recovered without any orders superpo- 
sition from the interference pattern. 

If we perform the hologram reconstruction by means of a plane reconstruction 
wavefront normal to the hologram plane and with the same amplitude of the plane 
reference beam, we obtain 


E(x, y) = 
= {ho + PT[E(x, y) + EDE, + Bree Es (x, y) + BTE e7 E (x,y) = 
= (ho + BrEZ)E, + PTE? (x, y)E, + BTEŽe P? Es (x, y) + BrEze "O € (x,y) 


The first term of the last expression is simply an attenuated version of the recon- 
struction beam and then, like this one, it propagates in the direction normal to the 
plane of the hologram, i.e. in the direction of the observer. The second term instead, 
while being proportional to the reconstruction beam, due to the spatial dependence 
of the coefficient E? (x, y), is composed of plane waves traveling at different angles 
with respect to the reconstruction beam propagation direction; however, if the fre- 
quencies of E?(x, y) are sufficiently small compared to the frequency of the carrier a, 
the energy of this term remains sufficiently close to the main propagation direction. 
The third term, the one giving rise to the real image of the object, is deflected by an 
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angle —@ with respect to the Z axis. Finally, the fourth term, the one giving rise to 
the object virtual image, is deflected by an angle +0 with respect to the Z axis. 
Thanks to this stratagem the beams are thus shifted each other and, for sufficiently 
high @ angles their overlapping is avoided, as shown in figure 1.3. 


Virtual image Real image 


Figure 1.3: Off axis holographic configuration. 


A completely analogous result can be obtained using spherical reference and re- 
construction beams with sufficiently large curvature radii and different curvature 
centres positions [5]. If the distance between the recording plane and the origin of 
the reference spherical wave is equal to the distance between the recording plane and 
the object position, we obtain a remarkable holographic configuration called lensless 
Fourier holography with peculiar properties. 


To estimate the value of the minimum angle at which we must operate in order 
to separate real image, virtual image and undiffracted term, it is sufficient to deter- 
mine the minimum value of the carrier wave frequency for which the various spec- 
tral terms do not overlap. In particular it is possible to demonstrate [5] that the three 
orders separation condition is fulfilled if the following inequality is satisfied 


8 > sin! 3BA 
Where B is the higher spatial frequency of the object wavefront. 


If the reference wave is much stronger than the object beam the zero order term can 
be ignored and the previous inequality reduces to [5] 


6 > sin"! BA 
In analog holography, where very high resolution recording devices are used and 


very narrow fringes are therefore recordable, it is possible to tilt the reference beam 
up to very high @ values and it is thus always possible to separate the three compo- 
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nents of the reconstructed wavefront. It is then possible to cut the undesired contri- 
butions by means of an appropriate spatial filter on the focal plane of a suitable con- 
verging lens placed in front of the transparency, that is on the plane where it materi- 
alizes the Fourier transform of the signal; after selecting the contribution of interest 
it is then possible to reconstruct the relative wavefront anti transforming the signal 
with a second identical lens. 


2. Digital holography 


From the point of view of the theory and the purposes, Digital Holography (DH) 
is directly derived from analog holography but differs from the classical technique 
for the hologram recording medium and for the wavefront reconstruction method. 
While, in fact, in classical holography the interferogram recording process is per- 
formed by means of appropriate photosensitive materials which need to be devel- 
oped before being used in the reconstruction process, in DH the recording medium 
is, more simply, a digital device, typically a CCD or a CMOS device. The interference 
pattern containing the information about the wavefront to be studied is thus sam- 
pled and digitized by the electronics of the device and stored in the memory of a 
computer. The reconstruction of the wavefront, obviously, can no longer be carried 
out in an optical way, at least not directly, but it is performed numerically by means 
of appropriate algorithms through standard calculation software; it is thus possible 
to derive the desired information on the wavefront under investigation, amplitude 
and phase, and to reconstruct the wavefront in a digital version too. 


2.1 Wavefront recording 


As we have just said, the recording process in DH is greatly simplified with re- 
spect to classical holography; analog photographic plates and complicated develop- 
ment techniques are no longer required but all it is needed is a detector working with 
the radiation in use and an image processing program to record the interference pat- 
tern. Current electronic recording devices, however, have a much smaller spatial res- 
olution with respect to old photographic plates: the resolution attainable with a clas- 
sic support can reach up to a maximum of about 7000 cycles/mm (where cycle 
means a pair of white/black fringes) while a CCD with square pixels of lateral di- 
mensions greater than 5 um (a typical values in CCD cameras) allows to solve a max- 
imum of 100 cycles/mm. This limitation is particularly restrictive when working in 
off-axis holography where the maximum recordable number of cycles/mm directly 
imposes a limit on the angle 0 between the reference beam and the normal to the 
sensor. If we refer, in particular, to speckle holography, in optimal recording condi- 
tions, the recording device should be able to solve entirely the interference pattern 
resulting from the superposition of the reference wave with all the waves spread by 
all the different points constituting the object under investigation. The maximum 
value of the angle @ is thus related to the lateral dimension of the object, to the lateral 
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dimension of the detector and to the distance between them. If, for simplicity, we 
refer to the simplified one-dimensional configuration sketched in figure 1.4, we can 
quantify these relations by means of simple geometrical considerations. 


Reference 
beam 


Object 


Figure 1.4: Maximum ð angle value according to the sampling condition. 


If we want that all the object point sources contribute to the interference pattern 
across the entire detector, it is necessary that the interference fringes between the 
reference beam and the highest spatial frequency component of the object wave can 
be recorded by the detector; if we indicate with @ the angle formed by the highest 
spatial frequency originating in the object with respect to the detector axis, we ob- 
serve [4], across the sensor, a sinusoidal interference fringe pattern between this 
highest spatial frequency and the reference beam with period 


À 


sin 0 + sin = 2 sin (=) eae (=) 


Which, for 6 = g and 0 + ọ = y, becomes 
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The maximum angle Ymax allowed by Whittaker-Shannon sampling theorem is 
reached when the value of the fringe period reaches its minimum admitted value 
Pmin equal to two times the detector pixels lateral dimension dp 


the maximum angle Ymax results thus to be 


a 
Wmax =2 sin! (=) 


For small y angle values, it is possible to approximate the sin function with its ar- 
gument, obtaining 


w small angle A 
max È 77. 

2d, 
Considered the typical pixel pitch values at each wavelenght, it follows, from the last 
equation, that the small angle approximation is usually well satisfied. 
With simple geometric considerations it is thus simple to calculate the maximum 
value max that is 
Dre? 
3; + 7 \ small angles A D+L 


a 
a = = sti (eee =t 283) Ee 
Onax = Ymax — P = 2 sin (È ) tan 2d, Fd 


p 


Where 

D is the object lateral dimension, 
L is the sensor lateral dimension, 
d is the object-sensor distance 


If the object lateral dimension is much larger than the lateral sensor dimension the 
previous expression simplifies further and becomes 


à D da-dyD 


max = 54 2d Zad 


p p 


This limitation on the angle 0, may lead, during the reconstruction process, to a 
partial overlapping of the three diffraction orders or, in other words, if the object 
dimensions are fixed, imposes a limit on the minimum distance between the object 
and the sensor: when the distance d is reduced the angle @ increases and, conse- 
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quently, the maximum admitted value of 6 decreases and may exceed the minimum 
value admitted to ensure the orders separation. 


Reference 
beam 


Object 


Figure 1.5: Minimum @ angle value to have the three orders separation. 


As we will see, in digital holography, the zero order term can be removed, during 
the reconstruction process, by means of various numerical algorithms and therefore, 
in order to separate the real and the virtual image, it is sufficient that the angle 0 sat- 
isfies the less restrictive inequality 


0 > Omin = sin! BA 


If we refer to the simplified configuration of figure 1.5, it follows that, for a fixed 
value of lateral object dimension and for a certain distance between the object and 
the recording device, the angle 0 can’t be smaller than the angle p between the high- 
est spatial frequency wave originating in the object and the normal to the sensor; in 
other terms, the minimum angle condition is obtained when the reference beam 
passes just adjacent to the object. 


2.2 Wavefront reconstruction 
As we have already mentioned, the wavefront reconstruction in digital holog- 


raphy is performed numerically by means of mathematical algorithms derived from 
the diffraction theory. 
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If we consider the hologram as an optical transmittance positioned on the re- 
cording plane (Xy,Yy), the reconstructed wavefront at a distance d, in the recon- 
struction plane (Xp, Yp), is in fact the result of the diffraction operated on the recon- 
struction beam by the hologram window (figure 1.6). 


Yo | Yu | YR 
b b b 
f Xo fn LAT f Xe 
ATC 
Optic axis Sa K | Ka 
| | — | de > 
Z 
y 
Object plane Hologram plane Reconstruction plane 
Z=-d Z=0 Z=d 


Figure 1.6: Object plane, hologram plane and reconstruction plane. 


To obtain numerically the analytical expression of such a wavefront we can ex- 
ploit the Rayleigh-Sommerfeld formula [7] that, by placing the obliquity factor equal 
to 1 [4], is written in this case 


PALI 
ize 


+00 
1 e 
E(Xp, Yr) = i Í R(x4,Yn)H (xy, Yy) dxydyy 


— 00 


p 


Where 

A is the wavelength used to reconstruct the hologram, 

H (xy, Yp) is the intensity of the interferogram in the hologram plane, 

R(Xy, Yy) is the wavefront reconstruction beam in the hologram plane, 

p = Jd? + (xr — xy)? + (Yr — yy)? is the distance between the generic point of the 
recording plane (Xp, Yg) and the generic point of the reconstruction plane (Xp, Yp), 

d is the distance between the recording plane (Xy, Yp) and the reconstruction plane 
(Xp, Yr). 


Once obtained the numerical expression of the complex reconstructed wavefront 
E(Xp, Yr), it is possible to extract from it the intensity value /(xg,yr) and the phase 
value $ (xp, Yr) of the object wavefront in the reconstruction plane 


mler) 


I(Xp, YR) = E(Xp Yr) E (Xe, YR) » r VR) = tan! ae 
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Where Re{£(xr,yr)} and Im{E (xp, yg)} denote, respectively, the real part and the 
imaginary part of the complex amplitude E (Xp, yp). 


The double integration of the Rayleigh-Sommerfeld formula involves considera- 
ble calculating effort but fortunately we will see that, in order to simplify its numeri- 
cal implementation, it is possible to convert this integral in one or more two- 
dimensional Fourier transforms; considered that, in case of digital recording, the 
complex function H (xy, Yp) in the hologram plane is discretized, the previous trans- 
forms are actually discrete Fourier transforms, more easily calculable by means of 
existing Fast Fourier Transform (FFT) algorithms. To operate this transformation 
there are two methods, the Fresnel method and the convolution method. 


Fresnel Method 
For values of Xp, Yr, Xy, Yy small compared to the distance d between the recon- 
struction plane and the hologram’ plane, the expression of 


p = Jd? + (xp — Xy)? + (Yr — yy)? can be approximated with its Taylor develop- 


ment 


Xp —Xy)? — yy)? Xp — Xp)? + — y,)2]2 
R a H) , Or Yu) _ [Gr H) Or- Yu)‘] Pa 
2d 2d 8d3 
The fourth term of this expression can be neglected if it is small compared to the 
wavelength [5], i.e. if 


[Œr - Xp)? + (Vr - Yu)? |? 3} [(xp — xy)? + Or- Yu)? 1? 
a” TE KA > d» Tata. ne 


Therefore, if we use the expression of the development up to the first order for the 
numerator (most critical factor) and up to the zero order term for the denominator 


(less critical factor) [5], we obtain, by developing the squares, 


E(xr, VR) = 


~ id vote) I Rx Yn)H (xy, Yn) etahi) eza C SHARE YD dyydyy 


This equation is called the Fresnel approximation of the Rayleigh-Sommerfeld inte- 
gral or Fresnel transform. If we now define the variables 


the previous integral becomes 
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E(u,v) = 
ia 90 
LIT 
E ame i R(Xy, Yy)H (xy, Yy) era ti) -inny dey dyy 
i 


The expression, in this way, has thus assumed, unless the multiplication factor out of 
the integral not depending on the variables xy,yy, the appearance of a two- 
dimensional Fourier transform and we can therefore write 


2 
ted 


ea R in 
E(u, v) = aa eitad(u*+v?)p {PR Cu: 4 )H Gen yy envio} 


Where F denotes the Fourier transform. 


It should be recalled, at this point, that, in digital holography, the hologram re- 
cording is performed in the digital domain and H (xp, yy) is therefore a discretized 
function; if we assume that the sensor is composed of a rectangular array of M by N 
pixels with spacing, along the axes XyandYy, equal respectively to Axy and Ayp, our 
hologram appears to be an array of numbers H (kAxy, l[Ayy) = H(k, 1) and the previ- 
ous integrals are therefore to be transformed into discrete summations or, equiva- 
lently, the continuous Fourier transform has to be replaced by a discrete Fourier 
transform; the wavefront on the reconstruction plane is then, in turn, a discrete 
function E(mAy, nAv) = E(m,n) of the discrete variables mAu, nAv; taking into ac- 
count that the maximum spatial frequency is determined by the sampling range in 
the spatial domain and namely [4] 


MAp=1/(Ax_H) , NA4v=1/(4y_H) 


It is possible to write 


Za 5 m2 n2 5 
oe aa Ri an zan DF {RCk, DH(k, peaaltari] 
U 


Where DF denotes the discrete Fourier transform. 


Finally, it can be observed that, according to Fourier transform relationship, 


Ad àd 


A = 
i YR NAYH 


and this means that the reconstructed wavefront in the plane (Xp, Yp), is represented 
by a matrix consisting of MXN elements, each of which is called reconstruction pix- 
el, with dimensions Axp, Ayp. It should be noted that the size of the reconstruction 
pixels depends on the reconstruction distance because the two expressions corre- 
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spond to the diffraction limit resolution, for, the optical system as a whole: the holo- 
gram constitutes the opening of this optical system and according to the theory of 
diffraction, the expressions of the values calculated above, Axz and Ayp, correspond 
to about half the diameter of the Airy disk of the diffraction pattern produced by this 
opening at a distance d [4]. 


If we suppose to have a square detector (M = N), with square pixels (Axy = Ayy = 
dy), we have 


a -24 
Pr Nd» 


Where Axr = Ayr = dpr 


From these expressions it follows that the reconstruction resolution increases 
when the distance d decreases; unfortunately this distance can’t be reduced indefi- 
nitely because there will be a minimum value d,,;, for which the maximum angle 
Omax admitted by the sampling theorem coincides with the minimum value of the 
angle Omin required to have the diffraction orders separation. In this condition, 
sketched in figure 1.7, we thus have the maximum resolution for a certain object di- 
mension. 


Reference 
beam 


Figure 1.7: Minimum distance configuration. 
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Working at the minimum distance we thus have 


Y 
8 = Omax = Omin = P = DT 
and, consequently 
a D+L È 2d,(D +L) D>L J 2d,D 
Ti — = — E 
4dp 2dmin ii ZI i ZI 


If we work at the minimum distance the reconstruction pixel pitch becomes 


2(D+L) 2D D»L 2D 
= =— + 2d — d,=— 


pr = N N p 


This means that, at the minimum distance allowed by the sampling theorem, the re- 
construction pixel pitch does not depend on the wavelength used to irradiate the ob- 
ject and can’t be smaller than the detector pixel pitch. Furthermore, if a large size ob- 
ject, with respect to the detector lateral dimension, is investigated, the reconstruction 
resolution increases with the number of elements constituting the detector. 


When a lensless Fourier holography configuration is employed the Fresnel 
transform simplifies significantly; the reconstructed wavefront in the sensor plane in 
this case is describe by 


+27 . 27 TT 
ewig Ve +%H?+yH° e ate aH tyn’) 
39) = -= E 


Vd? + xy?+yy? d 


When this expression is inserted in the reconstruction formula, Fresnel transform 
reduces to 


+00 


EQU, v) = Cemal498) fÈ Heey, yy) e OREM Ay dy = 
= Colma) P{H (xy, yy)} 
Where C is a complex constant. 
A lensless Fourier hologram is therefore reconstructed by a simple Fourier trans- 
form of the hologram function. Furthermore, when a lensless Fourier holography 


configuration is employed, it is no longer possible to focus on different areas within 
the object volume since the distance d does not enter any more the Fourier integral; 
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this means that in the reconstruction algorithm the distance d do not appear any- 
more and, consequently, in a lensless Fourier hologram reconstruction both diffract- 
ed orders appear in focus. 


Convolution method 
The convolution method consists in considering the Rayleigh-Sommerfeld for- 
mula for the wavefront reconstruction as a convolution integral, namely 


+00 


E(XR Vr) = | RH Vy) (xH, Yu) 9 XR — Xp Vr—Vu) dxudyy = [RH] * g 


=00 


Where 


1 pig VI tOR-XH)+OR-YH)? 
1) — Xy, YR-Yn) = = ee 
iñ [d2 + (xr — Xy)? + Or — Yu)? 
Applying the convolution theorem is then possible to write 

Ep yr) = F1{F{RH}F{9}} 

Similarly to the Fresnel method, to make the calculation in DH, we must write the 
discretized version of this expression and, if the analytical expression of the Fourier 
transform of the function g is denoted by G, we obtain 


E(m,n) = DF {DF {R(k, DH (k, D}G(m, n)} 


where [4] 


2 2 

M2 Ax p77 N2Ayy? 

2 H 2 YH 

A | p2 (mA (ta) 


G(m, n) = DF{9} = 3 a M2Axy? N2Ayy? 


This method does not impose a limit to the distance between the object and the 
sensor during the hologram recording process but requires the calculation of a Fou- 
rier transform and of an inverse Fourier transform. 

In this case the reconstruction pixels pitches, because of the double Fourier trans- 
form, are 


Axp =Axy , Ayr = Ayy 
These expressions could suggest that, by this method, the reconstruction pixels 


pitches have the same dimensions as the camera pixels, that is, for d > MAx?,/A, 
NAyj /A, smaller than the values obtained by the Fresnel method; actually the recon- 
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structed image resolution is still, at most, the resolution given in the formula of 
Fresnel method, as imposed by the diffraction limit, but, in this case, only a small 
portion of the object, equal to the size of the sensor, is played back; to obtain a com- 
plete reproduction of the object it is possible to use various techniques: it is possible 
to exploit the so-called zero padding technique, which substantially consists in add- 
ing a frame of zeros around the hologram matrix in order to perform the reconstruc- 
tion of a greater size matrix (the greater the synthetic matrix the greater the recon- 
structed portion of the object) [8;9]; alternatively it is possible to use a so-called nu- 
merical lens, i.e. to multiply the matrix of the hologram by an appropriate mathe- 
matical factor that reproduces the effect of demagnification of the field of view that 
would be achieved by the introduction of a lens in the real system [4;10]; another 
possibility is to reconstruct different portions of the object wavefront introducing a 
shift in the function G(m, n) and then stitching them together [4]. 


There exist different methods to remove the undiffracted term in DH [4]. As we 
have seen, the expression of the intensity recorded on the hologram is 


I(x,y) x Eo (x,y) + E? + 2E, (x, y)E, cosi, (x, y) — bo (x, y)] 


The first two terms of this expression give rise to the undiffracted disturbance 
term while the third term contains the information we are interested in. 

The most direct, but more laborious, technique to eliminate the effects of the 
first two terms of this expression from the reconstructed one, consists in acquiring 
three different images: the interferogram, the reference beam only and the object 
beam only; these images can be acquired from time to time by blocking the beam we 
are not interested in; subtracting from the image of the hologram the other two im- 
ages an undiffracted term free image is obtained. 

If we are not willing to acquire more than one image, we can use the subtraction 
of the mean value method: the third term of the reconstructed wavefront is statisti- 
cally varying between +2E,(x,y)E, from pixel to pixel at the recording device; the 
first two terms can then be suppressed by subtracting the average intensity of the 
hologram matrix from all the elements of the hologram matrix itself. 

A similar effect can be achieved if two holograms with stochastically changed 
speckle structures are subtracted from each other and the reconstruction of this dif- 
ference is performed. 

Alternatively, similarly to the case of analog holography, the elimination of the 
undiffracted terms can be obtained by cutting the undesired components from the 
spectrum of the signal; this result can be obtained in DH, numerically, in two ways: a 
first method consists in performing the Fourier transform of the signal, detect the 
noise components, remove them with appropriate filters and finally antitransform to 
obtain the hologram relative to the signal of interest. A second method consists in 
applying a high pass frequency filter or, preferably, a band-pass filter to the holo- 
gram itself, possibly after adjusting the parameters according to the information 
provided by the Fourier transform of the signal. The last method, based on more so- 
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phisticated filters, has been applied in all the hologram reconstructions performed in 
our experiments. 


2.3 Advantages of digital holography 


DH is characterized by a series of interesting characteristics that make it more 
appealing with respect to its classical counterpart and with respect to most standard 
imaging techniques; thanks to these characteristics it can be of interest in many dif- 
ferent applications and it is therefore meeting growing interest. 

The fundamental aspect that differentiates DH with respect to standard imaging 
techniques is the possibility to reconstruct both the amplitude and the phase images 
of the sample under examination; this information, of great interest in a wide variety 
of applications, is obtained with a minimally invasive technique since the remarkable 
sensitivity which characterizes modern recording digital devices allows, if necessary, 
to use very low power laser radiation during the hologram recording process. An- 
other interesting characteristic of DH is the possibility to perform, with a single ac- 
quisition, the reconstruction of the object image at different distances; this means 
that focusing on particulars belonging to different planes that make up the scene, 
can be obtained numerically during the reconstruction process starting from a single 
recorded hologram. Furthermore, as in classical holography, each portion of an hol- 
ogram carries the information of the entire sample and can be used to retrieve it; this 
means that even if part of the sensor is covered during the recording process it is still 
possible to recover the information about the entire scene under investigation. 

Many advantages characterize DH with respect to analog holography. First of all 
the digital hologram recording process is by far faster and easier than its analog 
counterpart: no delicate and time consuming photographic recording and develop- 
ment processes are required. Moreover, the times required for numerical reconstruc- 
tion are today, with a modern computer, in the order of a second and do are not a 
limitation in most applications. 

In addition to these important basic features of DH there is the possibility, 
through various techniques, instrumental or strictly numeric, to improve the quality 
of the final reconstructed image. It is possible, among other things, to improve the 
image resolution through the use of numerical zero padding technique [11] or 
through various different recording techniques [12], to enlarge portions of the image 
with the use of numerical lenses [10], to correct numerically aberrations of any lens 
used [13], to obtain an extended focus image of the sample [14], to correct the noise 
introduced by dust and inaccuracies in the alignment of the optics, by means of av- 
erage procedures over a large number of acquisitions. 
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Chapter 2 
Infrared digital holography 


The idea behind holography (both in its analog and digital version), that is to 
record and subsequently reconstruct the amplitude and phase information of the 
wavefront under investigation, is so general and rich in applications to be of interest 
in all regions of the electromagnetic spectrum. Obviously, the different wavelengths 
require different recording devices and the inhomogeneous technological develop- 
ment of the various acquisition devices has inevitably benefited certain regions of the 
spectrum with respect to others. In DH, for example, Infrared (IR) sources were pe- 
nalized compared to visible sources, due to the critical issues related to the detection 
of this kind of radiation. However, increasing interest toward standard IR imaging 
techniques, firstly related to their military, industrial and thermal efficiency applica- 
tions, has boosted the production of thermographic cameras and the availability of 
such IR cameras, based on the latest sensors generation, represents a significant 
stimulus to the development of DH techniques also in this region of the electromag- 
netic spectrum. This renewed interest could lead to innovative applications in secu- 
rity screening, night vision, and biological science, as it can extend holographic 3D 
imaging capabilities from visible to IR spectrum, up to the terahertz frequencies. 


1. Infrared radiation, sources and detectors 
1.1 Infrared radiation 


The IR range extends from about 0.750 um up to about 1000 um. The IR win- 
dow is usually conventionally divided into three different regions: Near IR (0.750 um 
- 3 um), Mid IR (3 um - 30 um), Far IR (30 um - 1000 um). These boundaries are in- 
deed not so firmly fixed and many other classifications do exist. The Mid IR region is 
sometimes called thermal Infrared (and also long-wavelength Infrared) while the Far 
IR region is nowadays commonly called Terahertz region (or T-ray region). Infrared 
radiation is used extensively in industrial and scientific applications such as envi- 
ronmental monitoring, industrial facility inspections, surgery, temperature sensing 
and spectroscopy. In most cases the above applications are based rather on Near and 
Mid IR radiation, due to the availability of efficient sources and detectors in these 
frequency ranges. On the contrary, the terahertz band, which lies between the optical 
and the microwave spectrum, is still almost unexploited and is currently one of main 
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research topics in optics because of its numerous and promising possible applica- 
tions. In this work we restricted the research field just to the longer wavelength of 
the Mid IR and Far IR radiation because of their attractive features with respect to 
DH employment. 


1.2 Sources 


The Mid IR experiments developed in this thesis were based on a well assessed 
Carbon Dioxide (CO2) Laser [15] because of the high versatility of this kind of coher- 
ent source: CO; lasers, in fact, are able to operate either pulsed or continuously with 
very high overall efficiency (up to 20%) and can produce continuous wave radiation 
ranging from few milliwatts up to more than 10 kW power (nowadays the highest 
power achievable with a commercial laser). They are based on a set of C0, vibration- 
al rotational transitions enabling them to emit from 9 um to 11 um; they can operate 
simultaneously on many lines or be tuned to a single wavelenght with very long co- 
herence lengths (up to several meters). The laser active medium is excited by an elec- 
tric discharge passing through the laser gas, which contains, besides carbon dioxide, 
also nitrogen (N2) and (usually) helium (He). Both N2 and C02 molecules absorb en- 
ergy from electrons in the discharge but the lowest vibrational level of N. easily 
transfers energy to C02. Helium helps to maintain the population inversion by facili- 
tating C02 molecules to drop from the lower laser levels to a lower level or to the 
ground state. Due to the high output powers that C02 lasers can reach, they are 
commonly used in cutting and welding industrial metal profiles and, due to the 
strong absorbance of water inside tissues at these wavelenghts, they are frequently 
employed in surgery and dermatology. A major drawback of this kind of lasers lies in 
their considerable size which, obviously, limits significantly the possibility to move 
out of the laboratory the techniques developed in this thesis using these sources. In 
this respect a promising alternative to CO; lasers is offered by a much more recent 
and extremely compact kind of lasers, operating both in the Mid IR and in the Far IR 
range, Quantum Cascade Lasers. 

The idea at the basis of Quantum Cascade Laser (QCL) working principle was 
firstly proposed in 1971 [16] but it was only in 1994, thanks also to the development 
of innovative atomic layer growing techniques, like molecular beam epitaxy (MBE) 
[17], that such a novel kind of semiconductor laser saw its first realization at Bell 
Labs [18]. Unlike typical interband semiconductor lasers which emit electromagnet- 
ic radiation through the recombination of electron-hole pairs across the material 
band gap, in QCLs the laser transition occurs between quantized electronic energy 
levels (subbands) within the conduction band of a repeated stack of multiple semi- 
conductor quantum wells; as a consequence, the emission wavelength is not deter- 
mined by a particular transition intrinsic to the material which constitutes the active 
medium, but can be tailored by acting on the quantum confinement characteristics; 
a broad wavelength range can therefore be covered, e.g., by tuning the thickness of 
the individual layers comprising the structure, without changing the materials com- 
position. QCLs are called unipolar devices in that only the electrons are involved in 
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the radiating process: every time an electron traverses the lasing medium it moves 
from one quantum well to the next, driven by the voltage applied across the device 
and, at precisely engineered locations, called active regions, it passes from one con- 
duction band energy state to a lower one thus emitting a photon. Additionally, while 
in semiconductor laser diodes, electrons and holes are annihilated after recombining 
across the band gap and can play no further part in photon generation, in a QCL, 
once an electron has undergone an intersubband transition and emitted a photon in 
one period of the structure, it can tunnel to the upper laser state of the following pe- 
riod where another photon can be emitted. This process of a single electron causing 
the emission of multiple photons as it traverses through the QCL structure gives rise 
to the term cascade. In this multistage cascade scheme, electrons are recycled from 
period to period, contributing each time to the gain and to the photon emission. 
Each stage belonging to the cascaded repetition consists of an active region followed 
by a relaxation-injection region. The electrons in the injection region are transferred 
by quantum mechanical tunneling to the excited state of an adjacent quantum well. 
Emission of radiation is a consequence of the electron relaxation to the ground state 
of the same quantum well (vertical transitions) or to that of the next one (diagonal 
transitions). To sustain population inversion, the ground level of the radiative transi- 
tion must be depleted by fast carrier transfer to the following electronic injector 
stage (i.e., via resonant optical phonon emission). Sequential repetition of these two 
regions allows electrons to be recycled from one stage to another, taking to the emis- 
sion, in principle, of as many photons per electron as the number of stages in the la- 
ser structure leading to a quantum efficiency much larger than one. 

Since the first realization in 1994, QCLs have experienced a remarkable devel- 
opment and, today, as previously stated, QCLs are fruitfully employed both in the 
Mid IR and in the Far IR, from 2.6 um [19] up to 250 um (1.2 THz) [20], with lasing 
demonstrated in both InGaAs / AlInAs and GaAs / AlGaAs material systems; an 
emission wavelenght encompassing a two orders range represents a really remarka- 
ble result, especially if it is considered that this broad operating range is obtained 
with the same lasing principle. It exists, however, a significant gap between Mid IR 
and Far IR QCLs performances. In particular, QCLs now represent the highest per- 
forming lasers available in the MIR frequency: they can operate almost continuously 
from 2.6 um to 24 um [21] delivering several watts of optical power in pulsed mode 
and several hundred milliwatts in continuous mode [22]; Mid IR QCLs can operate 
above room temperature in pulsed mode in a wide range of wavelenght (5 um - 12 
um) [23] and at room temperature in continuous wave for specific wavelenghts [24]. 

In the THz spectrum, QCL operation was firstly demonstrated in 2002 [25], and, 
at present, THz QCLs continuously cover the 1.2 - 4.5 THz range. Since their first 
experimental demonstration, THz QCLs have undergone significant progress with 
the development of specific active region designs (chirped superlattice, bound to con- 
tinuum and resonant phonon) but can still operate only at cryogenic temperature 
with a maximum output power of 250 mW, in pulsed mode, and 130 mW in contin- 
uous wave mode [26]. 
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1.3 Detectors 


As we have seen in the first chapter, array sensors are needed to perform DH 
and, among the most common IR sensors, we opted for pyroelectric and microbol- 
ometric sensors (most of which were already available in our laboratory). 

Pyroelectric sensors [27] are based on the ability of some materials (such as 
LiTaO3) to generate an electric potential when they are momentarily cooled or heat- 
ed; the temperature change induced in these crystals, slightly modifies the orienta- 
tion of the molecular electric dipole moments within the crystal structure, producing 
a variation of the polarization of the material; this variation, in turn, gives rise to a 
local electric potential that is possible to reveal and that disappears after a short re- 
laxation time; these sensors do not require cooling systems and are usable through- 
out the Infrared region. 

Microbolometric sensors [28] are based on the properties of some materials (such 
as amorphous silicon, aSi, and vanadium oxide, VaO) to change the value of its elec- 
trical resistivity by virtue of a rise in temperature; these sensors do not require cool- 
ing system and are able to work in a rather high range of wavelengths typically going 
from 8 um up to the Far IR with appropriate processing; thanks to their reduced di- 
mensions, microbolometric sensors allow to obtain very high resolution images, not 
achievable with other types of sensors and are therefore optimal devices for DH ap- 
plications. 


2. Advantages and disadvantages of IR DH 


Our decision to use long wavelenght radiation to perform DH is related to some 
intrinsic characteristics of this technique. First of all, as we have already mentioned, 
the radiation wavelenght does not enter into the expression for the reconstruction 
pixel pitch when operating at the minimum distance admitted by the sampling theo- 
rem; the reason for this unexpected property resides in the fact that, according to the 
general expression of the reconstruction resolution, the longer the wavelength in use 
the larger should be the reconstruction pixel size (the worse the resolution) but, if we 
look at the expression of the maximum angle between object beam and reference 
beam, we see that a longer wavelength allows to work at lower object-camera dis- 
tances thus compensating what was lost in terms of resolution. Reconstruction pixels 
pitch is indeed rather heavily dependent on the number and size of the elements 
constituting the recording device. In the past, Infrared detectors were characterized 
by a number of sensitive elements much lower than the detectors in the visible range 
and, at the same time, the size of the sensing elements was much larger than those of 
typical CCDs; probably this deficiency in IR detection devices, combined with the 
ease of use of visible radiation, prevented, in the past, a significant development of 
DH techniques at these wavelengths. Today the situation is greatly improved even if 
the gap between the performances of Infrared sensors and visible ones is not yet 
completely filled: Thermoteknix MicroCAM 1024HD, a high-level commercial un- 
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cooled long wavelenght detector, has a 1024x768 elements with 17 um of lateral di- 
mension, NEC IRV-T0831, a high level commercial uncooled Terahertz detector has 
320x240 elements with 23.5 um of lateral dimension, Hamamatsu ORCA-Flash 4.0, a 
high-level commercial visible camera, has 2048x2048 elements with 6.5 um of lateral 
dimension. With the expected progress of IR detection technologies, especially if we 
consider the significant current interest in the IR region, it is reasonable to assume 
that these values are going to be soon overcome so that the resolution gap between 
visible detectors and IR detectors are intended to get thinner and thinner. With these 
typical sensors technical values in mind it is simple to evaluate advantages and dis- 
advantages of Infrared DH (IRDH) with respect to visible DH: with reference to ta- 
ble 2.1 it is evident that, for a fixed object dimension D, IR radiation is the best 
choice if large size objects are of interest, even if a slightly lower resolution is ob- 
tained. 

As shown in the table, in fact, an object with 1m lateral dimension, for example, 
must be at least 24 m away from a visible detector while it could stay at about 3 m if 
CO2 radiation is used and even less than 0.5 m if we are working in the Terahertz 
region. The reconstruction resolution, instead, remains of the same order in all cases. 


Source CO: laser Quanta aes Solid State Laser 
Laser 
Wavelenght 10.6um 107um 0.532um 
Detector active pixels 1024x768 320x240 2048x2048 
Detector pixels size 17um 23.5um 6.5um 
Admin = (2d,D)/A 3.2D 0.42D 24.4D 
dpr = 2D/N 0.002Dx0.003D 0.006D x0.008D 0.001Dx0.001D 


Table 2.1: Typical hologram recording and reconstruction parameters in the Mid IR, in the 
Far IR and in the Visible region for a large size object of lateral dimension D. 


There is, however, also an intrinsic advantage in using long IR radiation in DH, 
and, more generally, in every interferometric experiment, that is the inherent lower 
sensitivity to vibrations of long wavelenght radiation. Optical path variations due to 
seismic noise or to any other vibration source, in fact, determine a lower phase varia- 
tion between the reference beam and the object beam, the longer is the wavelenght 
used to create the interference pattern. This benefit makes less critical the stability of 
the object and of the measuring apparatus in general, during the holograms acquisi- 
tion and is thus strongly indicated for DH, particularly in the investigation of large 
size objects. IR radiation makes it possible to work thus in less restrictive conditions 
allowing to operate even out of the laboratory, giving much more versatility to the 
technique and opening the way to many possible applications out of the laboratory. 
By means of long wavelenght radiation it is also possible to create video holograms 
(meaning holographic video of slowly changing dynamical scenes) without using a 
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pulsed laser and short acquisition times necessary if working with visible radiation. 
When very large size samples analysis is needed, CO; laser radiation is strongly rec- 
ommended since the high output power that CO: lasers can reach, is needed to irra- 
diate efficiently large object surfaces (such high power beams can in fact be expand- 
ed, through suitable optical systems, up to cover the entire surface of the object un- 
der study and ensure a sufficient power value all over the sample). Furthermore the 
very high coherence length typical of CO; lasers represents certainly an important 
advantage in every interferometric application. All these features make CO; lasers 
highly recommended if you need to use DH to investigate large size objects in out of 
laboratory applications. 

Among the reasons that may encourage the use of IR in DH there is also the 
transparency of various materials at certain wavelengths in the IR spectrum; this 
property can be exploited to study internal structures of such materials at different 
depths, to investigate heterogeneous materials or in security and safety applications. 
As we will see Mid IR radiation is highly transmitted through smoke while Far IR is 
well known for its capability of passing through plastic materials, clothes, paper, 
wood and many other materials. Unfortunately water is strongly absorbed both at 
10.6 um and in the Terahertz region and, consequently, exploitation of these wave- 
lenghts for biological purposes is restricted to thin samples investigation. 
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Chapter 3 
Mid IR digital holography experiments 


As anticipated, the first experimental part of this work was devoted to DH in the 
Mid IR and, more precisely, using the 10.6 um radiation produced by a CO); laser. 
Mid IR DH is gaining momentum thanks to its applications in the field of non- 
destructive testing analysis of large samples and in the field of real 3d imaging, and, 
as we will see, future applications could originate on the basis of some of the results 
obtained during this PhD project. 


1. Experimental apparatus 


As we are going to see in detail, the experimental apparatus needed to perform 
Mid IR DH is quite simple: all it is needed is a Mid IR coherent source, an appropri- 
ate array detector and few specific optics elements (beam splitters, lenses and mir- 
rors). Further instruments, like Spatial Light Modulators, motorized stages and spe- 
cial optics were employed for other applications and are described when needed. 

The laser used throughout our experiments was a CO; laser, Blade 100 by EL.EN 
Spa, emitting linearly polarized radiation at 10.59 um, relative to the line P(20), in 
continuous mode, on the Gaussian fundamental TEM00 mode, up to a maximum 
power of 130 W; the laser beam was characterized by a minimum diameter, beam 
waist, at the output of the resonant cavity, of about 10 mm and a divergence of about 
2.2 mrad. Using this kind of laser (belonging to the class IV) requires maximum at- 
tention both for the presence of a high power supply voltage (about 11 kV) and for 
the risk of eyes or skin exposure to direct beam or to scattered radiation by obstruc- 
tions on the path of the laser. To safely perform the alignment of the optical system 
we used a pointing source (laser diode at 670 nm wavelength) previously calibrated 
so as to follow the same path of the main laser beam. 

The detector used throughout our experiments was a microbolometric camera, 
by Thermoteknix, with the characteristics shown in table 3.1. The thermocamera, as 
in every holographic setup, was used in its lensless configuration. The software ac- 
companying this camera allows to acquire single images in different formats but also 
live video in the format *.avi. The acquired holograms, or video hologram, were nu- 
merically processed by means of ‘home made’ algorithms working in Matlab envi- 
ronment. 
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All tested configurations were arranged on a vibration isolated table by Melles 
Griot (used, for our purpose, not in vibration isolated condition). Since standard 
glass is not transparent at this wavelenght special Zinc Selenide (ZnSe) optics and 
gold coated plane mirrors were used. 


Model Miricle 307K — 25 
Detector material | Amorphous Silicon 
Array size 640x480 
Pixel pitch 25umx25um 
Spectral response 8um — 12um 
Sensitivity < 70mK 
Frame rate 50 frame/s 


Table 3.1: Thermocamera specifications. 


2. Experimental setups and results 


The very first results in IRDH at 10.6 um were obtained in 2003 by a INO group 
[29] using a CO» laser with 190 mW power and a pyroelectric camera, Pyrocam III 
by Spiricon, with 124 x 124 LiTaO; elements with pixel size 85 um x 85 um and cen- 
tre to centre spacing of 100 um x 100 um. In this first experiment a transmission 
holographic configuration was realized; with the setup illustrated in figure 3.1 the 
wavefront transmitted through a drilled metallic plate was recorded and numerically 
reconstructed both in amplitude and phase (see figure 3.2). 


CO? laser 


BSI [> Object beam 


Mirrer 2 


Figure 3.1: Transmission holography setup. BS1 and BS2 are 50/50 ZnSe beam splitters. 
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Figure 3.2: (From left to right) Hologram; amplitude reconstruction; phase reconstruction. 


In 2008 the same INO group, with a slightly more powerful CO2 laser (800 mW) 
and the same Pyroelectric camera, by means of a reflection holographic setup (figure 
3.3), managed to record and numerically reconstruct both the amplitude and phase 


of the wavefront reflected by small (1 inch diameter) metallic objects [30] (figure 
3.4). 


aw Ava 


Figure 3.3: Reflection holography setup. BS1 and BS2 are ZnSe beam splitters; M, M1 and M2 
are plane mirrors, BE is a ZnSe beam expander. 


Figure 3.4:.(From left to right) Amplitude reconstruction of an inscription on a metallic 
plate; phase reconstruction of a circular metallic disc with circular patterns inscribed. 
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2.1. First results 


When I started to work in this field, in collaboration with the same INO group 
which had realized the above experiments, our first efforts were dedicated to the im- 
provement of previously obtained results [31;32]; we thus decided to move to speck- 
le IRDH in order to investigate higher dimension scattering objects. As we have al- 
ready said, larger samples means longer object-detector distances and, consequently, 
lower reconstruction resolutions; in order to reduce this resolution loss we bought a 
higher resolution, larger size and more sensitive detector, the microbolometer detec- 
tor Miricle 307K described in the previous paragraph; in order to guarantee a more 
efficient object irradiation we opted for the 100 W CO; laser available in our labora- 
tory. The optical configuration used in these first experiments (figure 3.5) is a typical 
DH setup and was used as starting point for all the other configurations employed 
throughout our experiments; in what follows we will refer to this arrangement as to 
the basic configuration. 


-F 


CO, Laser | >f 


BS: M: 


`X 


Figure 3.5: Basic digital holography configuration. M1, M2 and M3 are plane mirrors; BS1 is 
a ZnSe 80/20 beam splitter; L1 and L2 are ZnSe 1.5 inch focal length spherical lenses. 


In this setup the CO; laser beam is first divided by a ZnSe beam splitter (BS1) 
which reflects 80% of the impinging radiation and transmits the remaining 20%. The 
transmitted part, which constitutes the reference beam, is reflected by means of a 
plane mirror (M1) towards a ZnSe variable attenuator (VA); the reference beam is 
then redirected by means of other two plane mirrors (M2, M3) toward the ther- 
mocamera but, before impinging on the detector, it encounters a ZnSe converging 
lens of 1.5 inch focal length (L2), which focus and then enlarges it in order to reach 
the thermocamera with enough low intensity and an almost planar wavefront. The 
reflected part of the fundamental beam constitutes indeed the object beam; before 
impinging on the sample it passes through a ZnSe converging lens of 1.5 inch focal 
length (L1), which focus and then enlarges the beam so as to irradiate a more or less 
large surface of the object depending on its distance from the sample. The interfer- 
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ence pattern created by the object beam and the reference beam across the sensor is 
optimized, in order to increase the fringes visibility, acting on the reference beam 
intensity by means of the variable attenuator. The hologram can be collected and 
digitally stored in a computer both in the form of a single image and in the form of a 
video if a dynamical scene is of interest. Thanks to the low sensitivity to vibration of 
the whole system the antivibration modality of the table is not activated during the 
experiments. Since the detector is only sensitive to IR radiation artificial lights or 
sunlight is not disruptive. Usually the sample is positioned just in front of the detec- 
tor and so great attention has to be done in the positioning of the mirror M3 because 
its position determines the angle between the reference beam and the object beam 
and, ultimately, the fringe spacing; as we have seen in chapter 1, it is necessary to 
vary the inclination between the two beams so as to obtain a sufficiently high fringes 
spatial frequency to separate the real image from the virtual one, without, however, 
breaking the limits imposed by the sampling theorem. When working with such 
high-power laser, special attention must be paid not to inadvertently send toward 
the sensor a too high reference beam which, considering the high sensitivity of this 
kind of detector, and its low damage threshold, would permanently injury it. 

In the first holograms recorded with this configuration we used, as object, a 
small bronze statue reproducing emperor Augustus (figure 3.6), about 10 cm high. 
In figure 3.6 the acquired hologram and its amplitude reconstruction are shown too. 
For completeness, the real reconstructed image, comprehensive of both diffraction 
orders, is shown in figure 3.7 (from now on we will omit to show both orders in next 
reconstructions). Both orders are in focus because, in this case, a lensless Fourier ho- 
lography configuration was used. 


Figure 3.6: (From left to right) Augustus small statue; Augustus IR hologram; Augustus nu- 
merically reconstructed hologram 


41 


Mid infrared digital holography and terahertz imaging 


Wavefront Amplitude Reconstruction 


100 200 300 400 500 600 


Figure 3.7:.Augustus numerically reconstructed hologram (both diffraction orders visible) 


A second and larger size statuette (34 cm high) reproducing Benvenuto Cellini’s 
Perseus sculpture has been successfully investigated with the same configuration; in 
figure 3.8 a photo of the statue, the relative hologram and the reconstructed ampli- 
tude wavefront are shown. 


Figure 3.8: (From left to right) Perseus statue; Perseus IR hologram; Perseus numerically re- 
constructed hologram. 


As we have anticipated, thanks to the lower vibration sensitivity of 10.6 um DH 
we were able to acquire, and successfully reconstruct, video holograms of dynamical 
scenes in CW laser condition. To test this possibility we positioned the sample on a 
motorized stage so that it could rotate on itself at low angular velocity (around 
360°/minute) and obtained a high quality video reconstruction (video 1). When 
higher angular velocities are employed the fringe visibility decreases progressively 
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and the reconstruction resolution gradually reduces; a higher frame rate and a short- 
er exposure time would be necessary to sample effectively the interference pattern 
and freeze the fringes in good visibility conditions. Our thermocamera could not ac- 
quire more than 50 frame/s but higher frame rate thermocameras (up to240 frame/s) 
are already commercially available. 


2.2 Super resolution by means of synthetic aperture 


Since our detector elements were larger and less numerous than typical visible 
detector elements our first results lacked of resolution with respect to typical visible 
holograms. One of the first improvements we decided to pursue consisted thus in 
increasing our holograms resolution [12]. Since hologram dimension (that is the de- 
tector surface) enters at the denominator in the reconstruction pixel pitch expression 
and since the detector surface is significantly smaller than the interferometric pat- 
tern surface created in air by the object beam and the reference beam interference, 
we decided to increase synthetically the numerical aperture of the system by means 
of an automatic technique capable of recording several shifted but partially superim- 
posable holograms and stitching them together. To achieve this result we exploited 
the basic DH configuration but fixed the thermocamera at two motorized transla- 
tional stages, LTS300 by Thorlabs (figure 3.9). 


M: 
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Figure 3.9:.Holographic setup for shifted holograms acquisition. 


By means of a remotely controlled routine it was thus possible to move the sensor 
horizontally and vertically along a serpentine path on the hologram plane and acquire 
various equally spaced, minimally overlapping, portions of the large interferometric 
pattern. The recording session had to be completed in a relatively short time so as to 
minimize changes in the interferential pattern even if the task was facilitated by the low 
vibration sensitivity of long wavelenght radiation. The numerous holograms thus rec- 
orded were then stitched together by means of an automated algorithm developed by 
our group, obtaining a digital hologram with a synthetic but larger numerical aperture. 
The stitching procedure, called registration, is a standard image processing technique 


43 


Mid infrared digital holography and terahertz imaging 


to determine the geometrical transformation that aligns points in one picture with cor- 
responding points in another picture. Often, the registration is performed manually or 
in a semi-automatic way by a user iteratively setting the parameters of the geometrical 
transformation. However, this approach is time consuming and can give subjective re- 
sults. On the other sides most of totally automated algorithms are not always applica- 
ble to align and stitch speckle holograms, which, due to their apparently random and 
very fine pattern and due to the inevitable inherently differences of the images struc- 
tures and tone dynamics, are hard to match and experimentally result in noisy images. 
We thus decided for a specific automated method based on the maximization of the 
mutual information (MMI) [33], which was capable of taking into account any possi- 
ble shift, rotation and scaling among the different holograms and that has shown ex- 
cellent results compared to the other methods. 

We decided to test this procedure to improve the reconstruction resolution of 
the previously investigated samples, Augustus and Perseus. In particular, a synthetic 
hologram composed of 7 x 3 standard holograms was employed to obtain the syn- 
thetic hologram and its relative amplitude reconstruction of Augustus image while a 
4 x 3 synthetic hologram was employed for Perseus image. The quite encouraging 
results are shown in figure 3.10 and 3.11. 

Again the high power of our CO; laser revealed as an advantageous parameter 
helping to obtain a larger interferential pattern with respect to standard visible pat- 
tern, and therefore allowing acquiring a larger number of single holograms to be 
stitched together. The reconstructed image resolution however cannot increase in- 
definitely because when the synthetic aperture hologram reaches a certain dimen- 
sion the reconstructed wavefront curvature, corresponding to different perspectives, 
cannot be properly represented with a 2D numerical reconstruction. The different 
perspective views written in the synthetic aperture hologram become evident if two 
different portions of the hologram (the left half and the right half part for example) 
are reconstructed and shown alternatively (video2). 


Figure 3.10:.Augustus hologram stitching. (From left to right) Single hologram and single 
hologram amplitude reconstruction 7 x 3; stitched holograms and stitched holograms ampli- 
tude reconstruction. 
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Figure 3.11: Perseus hologram stitching. (From left to right) Single hologram amplitude re- 
construction; 4 x 3 stitched holograms amplitude reconstruction. 


Noteworthy, the enhanced resolution makes it possible to exploit more efficiently 
the possibility, peculiar of DH, to focus on different object planes: thanks to the in- 
creased numerical aperture of the synthetic hologram, the depth of focus is decreased 
and it becomes possible to read the Augustus inscription on the base of the statue in 
the reconstructed image (left image of figure 3.12), or to focus on Augustus face (right 
image of figure 3.12), simply changing the numerical reconstruction distance. 


LI 


Figure 3.12: Stitched Augustus hologram amplitude reconstruction with Augustus inscrip- 
tion in focus (on the left) and with Augustus head in focus (on the right). 


2.3 Human size holograms 


As we have shown in chapter 2, by means of long wavelength DH, it is possible 
to obtain advantages in terms of larger field of view and lower seismic noise sensitiv- 
ity. Furthermore, using a high power source like our CO; laser, it is possible to irra- 
diate uniformly large surfaces. With long wavelenght radiation it becomes possible 
to record holograms in critical conditions including illuminated environments and 
not ground isolated setups in outdoor applications. Moreover, at this wavelength, 
holographic videos of moving objects, although with limited speed, can be easily 
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recorded and reconstructed. All these features are fundamental if large sample holo- 
grams are to be recorded. One of the peculiar difficulties in recording a large object 
digital hologram consists however in irradiating effectively the entire surface of the 
sample. In order to obtain optimal homogeneous sample irradiation we tested four 
different configurations and managed to record and reconstruct human size holo- 
grams [34] and, what is even more interesting, alive people holograms [35]. In order 
to compare the results obtained with each configuration, the object, a plastic manne- 
quin 1.90 cm high (figure 3.13), was maintained always in the same position with re- 
spect to the thermocamera. In particular, the distance between the object and the 
thermocamera was the minimum distance to record holograms of that size accord- 
ing to the equations discussed in chapter 1. The four setups shown below were just 
simple variations with respect to the basic configuration, differing from each other 
just in the way the sample was irradiated. 


Single beam setup 

The first configuration corresponded exactly to the base configuration except for 
the object beam lens Li which was substituted with a shorter (linch) focal length lens 
in order to obtain a larger beam size across the sample. This configuration however 
suffered from important issues: since the peripheral part of the gaussian object beam 
is characterized by a lower intensity, the object surface was not uniformly irradiated 
and this reflected in a poor and inhomogeneous reconstruction quality of the outer 
part of the sample. Further, the object beam dimension was limited by the focusing 
power of the lens and, consequently, only a portion of the sample could be fruitfully 
irradiated. In these conditions we could not take complete advantage of the in- 
creased field of view offered by the longer wavelength because the irradiated area 
was smaller than the maximum recordable size (figure 3.14). 


Infrared 
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Figure 3.13: The object, a plastic mannequin 1.90 cm high. Single beam setup. L1 is a ZnSe 1 
inch focal length spherical lens. 
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Figure 3.14: (From left to right) Mannequin image with irradiated area in red; mannequin 
hologram acquired with single beam setup; mannequin hologram amplitude reconstruction. 


Double beam setup 

A first solution to the inconveniences of the previous setup was obtained by 
means of the set up shown in figure 3.15. The object beam, after lens L1, was divided 
by means of a Zn Se 50/50 beam splitter, BS2, into two beams of equal intensity: the 
transmitted one was directed toward the lower portion of the sample, while the re- 
flected one was directed toward the upper part of the sample by means of the plane 
mirror M4. Using this configuration we were able to irradiate almost completely the 
mannequin with a more uniform energy density. As shown in the numerical recon- 
struction of figure 3.16, almost all the desired region was covered and the manne- 
quin became nearly completely visible in the reconstruction image. Similarly to the 
previous configuration, however, in those regions of the sample where the beams 
were less intense or superimposed, the reconstruction became less clear. 
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Figure 3.15: Double beam setup. M4 is a plane mirror, BS2 is 50/50 ZnSe beam splitter. 
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Figure 3.16: (From left to right) Mannequin image with irradiated area in red; mannequin 
hologram acquired with double beam setup; mannequin hologram amplitude reconstruction. 


Cylindrical lens setup 

In the third configuration (figure 3.17) a cylindrical lens was inserted just after 
the spherical one in order to obtain a stretched beam. This beam profile can be fruit- 
fully employed when the analysis of samples with a preferential direction, like our 
mannequin, is needed. Using the cylindrical lens it is possible to exploit more effi- 
ciently the available beam energy with evident advantage, as shown in figure 3.18. 
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Figure 3.17: Cylindrical lens setup. Levi is a ZnSe 5cm focal length cylindrical lens. 


Figure 3.18: (From left to right) Mannequin image with irradiated area in red; mannequin hol- 
ogram acquired with cylindrical lens setup; mannequin hologram amplitude reconstruction. 
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Scanning beam setup 

The last configuration (figure 3.19) is based on a completely different approach. In 
this case we used a single object beam, expanded by the usual spherical lens L1, and 
moved it along the sample, changing its propagation direction by means of the plane 
mirror M5 driven by two motion control devices (NewStep NSC200 by Newport). 


CO, Laser | = LA a 


Figure 3.19: Scanning beam setup: M4 and M5 are plane mirrors driven by two motion con- 
trollers. 


Figure 3.20: (From left to right) Mannequin image with irradiated area in red; mannequin 
hologram amplitude reconstruction at different scanning time and superposition of the most 
significant frames. 


Using this configuration we were able to uniformly irradiate, an overall surface 
of more than 4m’, an area larger than the admitted field of view. In this case we rec- 
orded a video hologram during which the object beam scanned entirely the manne- 
quin surface in about half a minute (video 3). When faster scanning acquisitions are 
realized the fringe visibility decreases progressively and the reconstruction resolu- 
tion gradually reduces; again this problem could be solved by means of a higher 
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frame rate and a shorter exposure time. Such kind of an acquisition procedure can 
be exploited in different ways: the video hologram can be reconstructed as a video 
which can be reproduced at normal velocity showing a slow scan of the sample, or it 
can be decimated and played at incremented speed so that the amplitude image of 
the whole object is perceived all at once; a third possibility consists in extracting the 
most significant frames and obtain the whole image by their superposition as shown 
in figure 3.20. 

The natural following step was to test the possibility to acquire alive human be- 
ing holograms [35]; using the first configuration we easily managed to acquire my 
half bust has shown in figure 3.21. Clearly, when operating with human beings, reso- 
lution get worse because of body micro movements and because of the lower reflec- 
tivity of human skin and clothes with respect to the plastic surface of the mannequin; 
it is, however, important to point out that such a result would be very difficult, if not 
impossible, to reach with a CW visible laser in not isolated conditions. 


Figure 3.21: Amplitude reconstructions of my half bust with single beam setup. 


Holography and, more generally, interferometry are commonly used for non- 
destructive testing of mechanical structures, stress analysis and quality control. The 
advent of DH has represented an invaluable advance in this field since it offered the 
possibility to perform real time investigations with a less expensive technique. A 
strong limit to this kind of application, especially for industrial purposes, is however 
represented by the presence of vibrations and by the limited sample dimensions ad- 
mitted in visible DH. As we have demonstrated, IRDH, thanks to its capability to in- 
vestigate large size samples without any particular attention to vibrations, represents 
thus a perfect candidate to substitute visible radiation for such purposes. Recently, 
also in the wake of the results obtained with our works, other research groups have 
developed such kind of non-destructive analysis by means of Mid IRDH [36]. With 
this technique it is in fact possible to follow, in real time, large samples deformations 
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and quantify their amount evaluating the induced phase variation. In figure 3.22 a 
thermal expansion (due to laser heating itself) of Perseus statue has been revealed as 
phase variation with respect to the undisturbed initial sample condition (video 4). 


Figure 3.22: Reconstructed phase changes during Perseus statue heating with respect to the 
undisturbed initial condition. 


2.4 Visible analog reconstruction of IR digital holograms 


Holograms reconstruction in DH is usually implemented by means of numerical 
algorithms but nothing prevents to ‘write’, somehow, the hologram on an appropri- 
ate medium in order to perform an analog reconstruction. Spatial light modulators 
(SLM) (arrays of pixels where each pixel modulates the phase and the amplitude of 
the light transmitted through or reflected from it) are the elective devices to accom- 
plish this particular task [37;38]. However, in general, the SLM pixel pitch dpsım dif- 
fers from the digital hologram (recording device) pixel pitch d, and this discrepancy 
may represent an important complication. Furthermore an optical reconstruction 
makes sense only if visible light is used while our holograms were recorded in the IR 
range. The so called imaging equations [4;5] relate the coordinate of an object point 
with that of the corresponding point in the reconstructed holographic image, in the 
most general conditions, that is for a rescaled hologram and for different record- 
ing/reconstruction wavelengths. According to these equations, the distance at which 
the reconstructed image appears in focus is [4;5] 


1 uf1 LA: 
Zim = | —+-5(—- 
Zrec M° \Zref Zobj 
Where 


Zim is the Z coordinate of the reconstructed object image, 
Zop; is the object distance from the hologram plane, 
Zref and Zrec are the coordinates along the Z axis of the origin points of the spherical 


reference and reconstruction beams, respectively, 
dpSLM . 
= a is the scale change factor of the hologram, 
p 


51 


Mid infrared digital holography and terahertz imaging 
u= “2 and A, and A, are the recording and reconstruction wavelengths, respectively. 
sl 


The upper set of signs applies for the one image wave and the lower set for the other. 
When Zim is negative, the image is virtual and lies to the left of the hologram, while 
when it is positive, the image is real and lies to the right of the hologram. Depending 
on the geometry it is possible for one image to be real and for the other to be virtual 
or for both to be real or both virtual. 


Furthermore, according to the imaging equations, the reconstructed image under- 
goes a transverse magnification and longitudinal (or axial) magnification respective- 
ly equal to [4;5] 

2 


M,=m(1¥2 sob) ont” M, = =M? 


H Zrec Zref pei u 
We note that, if the hologram is recorded in a lensless Fourier holography configura- 


tion (Zon; = Zref )» 


Zim = Zrec 


regardless of the wavelengths ratio and of the scale change factor of the hologram 
and the magnification factors simplify as follows 


2 
tT Zi Zi 
M,=F4(2m), m= (1) 
m Zobj m Zobj 


Generally, therefore, when a 3d hologram reconstruction with a different wave- 
length is performed, the resulting image is affected by aberrations (spherical aberra- 
tion, coma, astigmatism, field curvature and distortion). The expression of these var- 
ious aberrations can be deduced from the imaging equations and it can be demon- 
strated [39] that lensless Fourier holography configuration minimizes most of them. 


Single SLM 

In a first work [40], 10.6 um holograms of the habitual Perseus statuette were 
recorded in the basic holographic arrangement but in a lensless Fourier holography 
configuration. In particular, we acquired 120 holograms rotating, each time, the fig- 
urine around itself of 3 degrees. In the reconstruction process, we used a diode 
pumped solid state laser emitting at 0.532 um. The laser beam was optically manipu- 
lated in such a way as to obtain a converging beam impinging on a reflective liquid 
crystal on silicon (LCoS) phase only SLM (PLUTO by Holoeye with 1920 x 1080 
pixels, 8 um pixel pitch and 60 Hz frame rate) where the acquired holograms were 
played back in sequence with regular frequency (figure 3.23). 
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Figure 3.23: Experimental setup for recording-reconstruction process. M is a mirror, L is a 
lens, BS is a beam splitter, MO is a microscope objective, SLM is a spatial light modulator, SF 
is a spatial filter. 


With such an arrangement it was possible to obtain an analog reconstructed video 
of the rotating figurine on a screen at a fixed distance from the SLM (video 5). An im- 
age of the optically reconstructed wavefront, acquired by means of a standard CCD 
camera, is compared with a standard numerical reconstruction image in figure 3.24. 


Figure 3.24: (a) Numerical reconstruction of Perseus statuette hologram (b) SLM optical re- 
construction of Perseus statuette hologram performed at a visible wavelength and collected by 
a monochromatic CCD. 
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The experimental reconstructed image position and its magnification with re- 
spect to the real object, were in good agreement with the values predicted by the the- 
oretical expressions and no significant aberration was observed. In order to verify 
that no aberration really affected the displayed image in our case, we inserted our 
experimental conditions values in the theoretical aberration coefficients expressions 
confirming that the wavefront aberrations were all negligible (<1 wavelength) and 
did not substantially affect the reconstructed images. 


Finally, by means of a numerical procedure [41], we created synthetic scenes that 
could be optically displayed. The procedure is based on recording several digital hol- 
ograms of individual objects covering a 360° perspective. Each hologram is geomet- 
rically transformed to change the object’s position and size within a very large depth 
of field and it is thus possible to create a dynamic 3D scene in which the object trav- 
els backwards and forwards in the 3D volume while performing a “pirouette”. Com- 
bining holograms corresponding to different objects, it is also possible to create dy- 
namic scenes with two or more objects rotating while travelling back and forth (fig- 
ure 3.25). 


Figure 3.25: Numerical (a, b) and optical (c) reconstructions of the synthetic holograms ob- 
tained by combining the digital holograms of Perseus and the digital holograms of a similar 
statue representing Botticelli’s Venus. The reconstruction distance is chosen in order to have 
Perseus and Venus in focus when they are close-up (a) or when they are far away (b). 


Multiple SLMs 

In a second work [42], in collaboration with the Department of Electrical and 
Electronics Engineering of Bilkent University (Turkey) and with the Bulgarian 
Academy of Sciences of Sofia (Bulgaria), we managed to obtain a holographic digital 
video display of a three-dimensional ghostlike image of Perseus multiview recorded 
holograms, floating in space. To achieve this result we used a holographic video dis- 
play system built from nine LCoS phase only SLMs (Holoeye HEO-1080P with 1920 
x 1080 pixels, 8 um pixel pitch and 60 Hz frame rate) forming a circular configura- 
tion (figure 3.26). 
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Figure 3.26: Circular holographic display. Left, arrangement of the nine phase-only SLMs, 
denoted as 1...9. Right, illumination of a single SLM. 


Elimination of the gaps between the SLMs was provided by a beam splitter to tile 
them side by side [43] and to achieve a virtual alignment with a continuous in- 
creased field of view. To position the reconstructed 3D image slightly above the dis- 
play setup and to avoid blocking of the observer’s vision by the display components, 
the SLMs were also tilted up at a small angle. Negligible reduction in the quality of 
the reconstructions for a tilted illumination of up to 20° has been shown by experi- 
ments [44]. All SLMs were illuminated with a single astigmatic expanding wave by 
means of a cone mirror as shown in figure 3.26. With the help of this configuration 
observers can see 3D ghost-like image floating in space and can move and rotate 
around it (figure 3.27 and video 6). 


Figure 3.27: Left, single SLM optical reconstruction at 0.532 um of the hologram captured at 
10.6 um and projected on a screen. Right, ghostlike multiple SLMs optical reconstruction of 
the same hologram. 


These two works demonstrates the possibility of obtaining direct real time 3D vi- 
sion of IR recorded digital holograms and, considered the previously underlined ca- 
pability of IRDH to investigate large size samples, emphasize the potentialities of 
IRDH as valuable candidate in the research field of real 3d TV and for possible virtu- 
al museum applications. 
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2.5 Smoke and flames hidden object holograms 


In this section we report our latest results in the field of IRDH at 10.6 um; in par- 
ticular we illustrate one more useful peculiarity of this technique, that is the possibil- 
ity to detect object, or even animate beings, through smoke and, what is even more 
noticeable, through a wall of flames. The ability to monitor objects completely hid- 
den, in the visible range, by smoke and flames and the even more interesting possi- 
bility to detect alive people in fire scenarios is a key and challenging target for its 
possible industrial applications and for its possible remarkable military and, above 
all, safety purposes. Visible radiation is strongly affected by smoke or fog and vision 
can be completely impaired in smoky environment. On the contrary last generation 
of uncooled IR microbolometer detectors, commercially available for imaging in the 
thermal Infrared, allow passive or active (i.e. with IR laser illumination) clear vision 
through such scatterers since the IR electromagnetic radiation is just slightly scat- 
tered by fog drops and smoke particles. This property of IR radiation is well known 
and in fact many fire departments already use this kind of technologies for exploring 
fire scenarios in order to have a clearer vision and to operate safely in such hostile 
environments. Unfortunately, in presence of flames, even these detectors cannot be 
of great help since the electromagnetic radiation emitted by flames can severely satu- 
rate them, as well as standard CCD or CMOS, occluding the scene behind the 
flames. 

As we have repeatedly stated, DH at long IR wavelength has considerable ad- 
vantages, with respect to visible DH, that make it very flexible and useful for record- 
ing real world large scenes. Thanks to the longer wavelenght, IRDH is characterized 
by a lower seismic noise sensitivity and a larger field of view; as we have proved, the- 
se features allow to record digital holograms of large dynamical scenes and make it 
possible to obtain human moving holograms in CW conditions; the broad availabil- 
ity of high power CO2 laser sources with good coherence properties makes it possi- 
ble to expand the object beam and achieve an efficient irradiation across very large 
scenes without too much restrictive limits on the optical path difference between ob- 
ject beam and reference beam. Furthermore IR digital holograms can be recorded in 
open space and daylight conditions thus allowing this technology to be directly 
brought out of the laboratory. Another important characteristic to remember is the 
possibility to use only a portion of the acquired hologram to reconstruct the entire 
scene even if at the expenses of resolution. Finally DH usually employs a lensless set- 
up to obtain out of focus acquisitions and to recover numerically, in the reconstruc- 
tion step, the object wavefield at the desired focus plane. 

All these features, combined with the possibility to exploit IRDH to see behind 
smoke and flames, which is, to our knowledge, not achievable by any other imaging 
technique, provides the unique possibility to perform real time dynamic detection of 
moving people in fire scenes. Indeed, as we have proved, IRDH allows not only to 
record objects immersed into smoke but also to visualize objects behind flames in- 
dependently from the chemical nature of the involved burning materials and from 
their emission spectrum. On the contrary, we have verified how standard IR imaging 
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cameras cannot carry out this task. In order to demonstrate these unique imaging 
potentialities, two kinds of experiments were conducted, where the capability of see- 
ing through smoke and flames by IRDH were respectively tested. 


Seeing through smoke 

In a first series of experiments [35] we tested the capability of IRDH to work ef- 
ficiently in smoky environments. In order to investigate this condition we used the 
basic holographic setup but immersed the object under investigation in a thick blan- 
ket of smoke (figure 3.28). As test object we used the habitual small Augustus statu- 
ette but, in order to obtain enough high smoke density around it, we put it inside a 
sealed Polymethyl methacrylate (PMMA) box. In one face of the box two windows 
in the Infrared range where fixed: an input AR/AR ZnSe input window through 
which the laser beam could reach the object and a germanium output window 
through which the light diffused by the object could reach the thermocamera. 
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Figure 3.28: Imaging through smoke setup. 


By means of a lateral aperture practiced in the box we let the smoke obtained by 
burning incense in a small furnace, enter into the box and, at the same time, we 
started acquiring a holographic video of the interferometric pattern with the increas- 
ing smoke density and a standard video in the visible range to testify the complete 
invisibility of the sample at the maximum smoke density (video 7). In order to have 
a quantitative measurement of the smoke concentration inside the box at any mo- 
ment, we evaluated the intensity attenuation ofa 15 mW laser diode radiation travel- 
ling 6 cm inside the box before reaching a silicon photodiode; by means of an HP 
3401A multimeter we could check the current at the photodiode output so that the 
extinction coefficient at any smoke density could be evaluated by means of Lambert 
Beer law. 

In figure 3.29 the visible images of the empty box and after inletting smoke are 
shown; evidently at high smoke density the vision is completely impaired because of 
the severe scattering of the visible radiation; at the maximum smoke density the cur- 
rent in the photodiode was reduced by two orders of magnitude and it was impossi- 
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ble to see the statuette even from the short side of the box. In figure 3.29 the images 
obtained with the thermocamera working in normal mode with its IR objective in 
the two conditions are also shown. As expected, with a standard thermographic im- 
age the object immersed in the smoke is still clearly visible since IR light at such long 
wavelength is only slightly scattered by smoke particles (according to Rayleigh and 
Mie scattering laws). Finally, the reconstructed IR holographic images in the two 
conditions, shown in figure 3.29, demonstrate that a clear vision through smoke can 
be obtained by means of IRDH too. 


Figure 3.29: (a) Augustus statuette in the PMMA box without smoke. (b) Augustus statuette 
in the PMMA box with smoke. (c) PMMA box without smoke: thermographic imaging of 
Augustus on the left and hologram amplitude reconstruction on the right. (d) PMMA box 
with smoke: thermographic imaging of Augustus on the left and hologram amplitude recon- 
struction on the right. 


It is important to note here that the random movements of the scattering particles 
(smoke and dust particles), which represents a noise source in the standard thermo- 
graphic image, on the contrary may contribute to reinforce the clear vision in the hol- 
ographic image: indeed we will show that, if multiple acquisitions are reconstructed 
and opportunely averaged, improved images with higher resolution can be obtained. 


Seeing through flames 

In a second series of experiments [35] we tested the capability of IRDH to detect 
an object hidden behind a curtain of flames. In order to investigate this possibility 
we used the basic holographic setup but inserted flames between the thermocamera 
and the object under investigation (figure 3.30). 
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Figure 3.30: Imaging through flames setup. 


The first tests were performed using as object the habitual Augustus statuette but 
then the setup was rearranged in order to test the technique with human beings. The 
flames were obtained employing candles in case of Augustus holograms and portable 
mini stoves for human size targets. In both cases a large portion of the object was 
covered by flames, impairing its vision to a normal CCD camera but also to the 
thermocamera. On the contrary, the holographic recorded images showed no such a 
problem, allowing to see through flames without significant resolution loss. 

This capability can be explained by the fundamental intrinsic features of holog- 
raphy. First of all, since no objective is required during holograms recording, the IR 
radiation energy emitted by the flames is not focused on the detector but is distribut- 
ed over its whole surface; this means that no image of the flames is formed on the 
detector and, consequently, no pixel saturation effect is obtained (figure 3.31). 


Figure 3.31: (Top) Detector saturation in presence of the thermocamera lens. (Bottom) Ho- 
mogeneous background in the lensless holographic configuration. 
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In other words, since out of focus images are recorded in DH, the typical satura- 
tion effect of the IR camera elements observed in standard imaging configurations is 
avoided and the sensor is not blinded by the flame emission. If the radiation energy 
intercepted by the sensor should be so high to disturb, somehow, the detector, it 
would be however always possible to use a narrow band filter around 10.6 um to 
completely remove its contribute. 

A second important aspect of the question, strictly related to the interferometric 
character of DH, is that the flames radiation is not coherent with the radiation used 
to get the interferogram and consequently it does not affect the interferometric pat- 
tern in any way. 

Furthermore, thanks to the ability of holography to reconstruct entirely the ob- 
ject wavefront from a smaller portion of the hologram, this technique allows an en- 
hanced vision even if some macro-particles, often present in real fire-scenarios, ob- 
struct direct imaging. 

The only difficulty in recovering the object image with this technique could be 
represented by the refractive index changes due to flames gases vortices which may 
affect the object wavefront in the same way in which hot air vortices in front of an 
object would distort its image in the visible range; the complication could originate 
from an excessively fast changing rate of the refractive index and, consequently, of 
the wavefront, but, even in this case, this complication could be solved with a higher 
frame rate thermocamera and lower exposure times. 

The advantages of digital holographic image system with respect to standard 
thermographic and visible recording systems are clearly appreciable from the com- 
parisons of the visible, the thermal and the holographic reconstructed images of the 
samples in figures 3.32, 3.33 and 3.34 (and video 8). 


Figure 3.32: Imaging of Augustus statuette seen through candles flames. (Left) Thermo- 
graphic image. (Right) Holographic reconstructed image. 
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Figure 3.33: Imaging of a hand seen through stove flames. From left to right: visible image, 
holographic reconstructed image, enhanced holographic reconstructed image (see the end of 
the chapter for details on the enhancement procedure). 
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Figure 3.34: Imaging of an alive human being seen through two stoves flames. From left to 
right, thermographic image without flames, thermographic image with flames, visible image 
with flames, holographic reconstructed image with flames. 


We hope that these preliminary but promising experimental results could repre- 
sent an important boost toward a detailed experimental analysis of the applicability 
of this technique in the industry and safety fields. 


2.5.3 Speckle contrast reduction by multiple acquisitions 


A further advantage of the proposed technique is that, through a numerical pro- 
cessing of the acquired holograms, it is possible to improve a posteriori the image 
quality reducing speckle noise. This feature could be exploited if, for example, it 
should be necessary an a posteriori analysis of the fire scene. 

Speckle is a multiplicative noise which occurs whenever a coherent light hits on 
targets whose roughness varies on the same scale of the employed wavelength. In 
this case, each detector element records the coherent superposition of a number of 
different scattering contributions. As they experience microscopically different 
paths, their phases can be very variable and the result of the coherent superposition 
to the receiver is the typical succession of dark and bright spots known as speckle 
noise. As a result, the quality of the reconstructed hologram in DH gets worse in 
terms of contrast and pixel resolution. Due to the nature of the problem, it has been 
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tackled in literature relying on statistical approaches. A common way to get indica- 
tions about the speckle extent is to measure the speckle contrast as follows: 


lol 
C== 
u 

where o and u respectively denote the standard deviation and the mean intensity 
of the image. 


Furthermore, a punctual measurement of the intensity variations due to the 
speckle noise is obtainable by calculating the relative deviation as 


I(x,y)-I 
Roy) = E2 


Where /(x, y) is the intensity of the pixel (x,y) of the reconstructed hologram 
and / its mean intensity calculated over the whole image. In particular, if a homoge- 
neous portion of the original image is chosen, a smooth behaviour of the plot should 
be expected. So, any rapid variations and sudden spike in its shape have to be at- 
tributed to the speckle. 

In order to reduce the speckle contrast it is possible to exploit some kind of di- 
versity during holograms recording to obtain holograms from a set of uncorrelated 
speckle patterns and average them to get a multi-look (ML) image. Indeed, it can be 
shown [35] that the superposition of N statistically independent random variables 
returns a random variable whose variance is reduced of a factor 1/N . Hence, the 
speckle contrast improves of a factor 1/VN with respect to the contrast of a single- 
look (SL) image. 

Figure 3.35 shows the holographic reconstruction of Augustus statuette placed 
into the PMMA box of figure 3.29. In particular, we acquired a set of holograms after 
inletting smoke into the box, whose effect is to provide a temporal diversity. After a 
proper decimation of the original hologram stack we got a subset of N = 25 holo- 
grams and for each of them we performed the numerical propagation after dropping 
the diffraction orders out of interest. The reconstructed amplitudes have been aver- 
aged and speckle measurements have been carried out on both the SL and the ML 
images in terms of the contrast estimators defined above. 

Figure 3.35 (a) shows the SL image and the relative deviation calculated over the 
homogeneous area corresponding to the red box in figure. As expected, the noise is 
responsible for the rapid fluctuations of R(x, y) and the degradation due to speckle 
noise is evident in the amplitude image as well. In figure 3.35 (b) the described pro- 
cedure is sketched and the resulting ML output is shown. The ML improvement is 
clearly appreciable in the amplitude image and the relative deviation exhibits a much 
smoother behaviour. This suggests that a significant gain has been achieved in re- 
ducing the speckle noise by combining multiple acquisitions. This claim is con- 
firmed by the measure of speckle contrast, which gets reduced of the 12% with re- 
spect to the SL image. However, a residual correlation between the employed frames 
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is present, which in turn results in the gap with respect to the maximum theoretical 
improvement factor, i.e. the 20% in case of N = 25. As a further validation, we cal- 
culated the amplitude histograms in a homogeneous test area of both the SL and the 
ML reconstructions. As expected in the ML gain traduces in a smaller variance, 
which is apparent in the plots shown in figure 3.35 (c). 


Figure 3.35: Imaging through smoke. (a) Holographic amplitude reconstruction before nu- 
merical processing and Relative Deviation corresponding to a homogeneous portion of the 
image (red box in figure). (b) Speckle reduction by processing a time sequence of holograms: 
multi-look reconstruction and Relative Deviation improvement. (c) Amplitude histograms: 
comparison between the single look (upper curve) and the multi-look image (lower curve). 


In case the line of sight between the target and the recording device is impaired 
by the presence of a flame, the gas vortices provide in turn a temporal diversity and a 
set of uncorrelated holograms can be acquired and temporally averaged to get again 
an improved ML output. Figure 3.36 shows the improvement achievable by numeri- 
cal processing. In particular, in the SL image the effect of the saturation is apparent, 
with dark areas hindering a clear vision of the human target behind the flames; the 
human target becomes clearly visible in the ML reconstruction. In this way more de- 
tails can be appreciated and a post analysis of the room scene is possible. 
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Figure 3.36: Imaging of a human target behind a flame. Left, SL holographic reconstruction. 
Right, ML holographic amplitude image. 
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Chapter 4 
Terahertz imaging experiments 


The Terahertz frequency range has historically been characterized by a relative 
lack of convenient radiation sources, detectors and transmission technology and, for 
this reason, it has remained one of the least developed spectral regions. Over the last 
few years, however, THz radiation has become of large interest for its peculiar sensi- 
tivity to many molecular absorption lines, appealing for spectroscopic, astrophysics 
and biological analysis, and for its potential in non-invasive imaging applications, 
spanning from medical diagnostics to homeland security. This renewed interest has 
benefited particularly of the invention of Quantum Cascade lasers (QCLs) in 1994 
[18] as well as of the recent development of electronic emitters, such as resonant 
tunnelling diodes, and frequency multipliers [45]. In comparison with visible or 
Near Infrared waves, THz radiation can penetrate organic materials such as plastics, 
clothes, wood, paper products and thin biological samples. Similarly to X rays, THz 
radiation cannot penetrate metals but differently from them, it is a not-ionizing ra- 
diation and, thanks to its low photon energy, it does not cause damage to the investi- 
gated materials. Thanks to such peculiar properties, THz radiation can thus be fruit- 
fully employed for security screening and illicit material detection (safe body scan- 
ners, packaging inspection), non-destructive testing and quality control. Since the 
onset of THz imaging, in 1995, by Hu and Nuss [46], who coined the term T-rays, 
various imaging modalities for numerous applications have been proposed, devel- 
oped and demonstrated. At the beginning the majority of research systems in THz 
imaging made use of Time Domain Spectroscopy (TDS) but several other configura- 
tions, including time of flight imaging [47], tomographic imaging [48], real time ar- 
ray detectors imaging [49;50], near-field sub wavelenght resolution imaging [51;52], 
have been successfully applied so far. 

The second part of this PhD project has been devoted to this continuously de- 
veloping and promising field of THz imaging. In particular we developed a simple 
transmission imaging configuration which worked as starting point toward the al- 
most unexplored field of THz digital holography. 


1. Experimental apparatus 


In all the experiments performed in the Far IR we employed the coherent radia- 
tion emitted by a Quantum Cascade Laser (QCL) processed at the Nanoscience In- 
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stitute of Pisa (CNR); in particular the source at our disposal was a bound to contin- 
uum [53] QCL emitting at 2.8 THz (107 um) driven in pulsed regime. The laser 
could operate up to 85 K and, consequently, a cryostat chamber (Ricor LV FA Cool) 
and a vacuum pump (Pfeiffer Vacuum Hicube TC110), to create vacuum in the cry- 
ostat chamber, were employed. A waveform generator (Agilent 33250A) combined 
with a current driver (Wavelength Electronics QCL 1000, 0.2 A/V transfer function) 
were used to drive the laser with the desired pulse duration. In graph 4.1 the QCL 
power vs current characteristic (5 Hz square wave driving current with 50% duty cy- 
cle at T = 20 K) is illustrated. As shown in the graph, in this operating conditions, 
the laser threshold is reached at about 0.6 A and the maximum laser power (about 
1.5 mW) is obtained at about 0.9 A of driving current. 
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Graph 4.1: Laser output power vs laser driving current. 


To detect laser radiation we had available in the laboratory different detectors: 

- Pyroelectric array detector Pyrocam III by Spiricon (124 x 124 LiTaO; ele- 
ments, pixel size 85 um x 85 um, centre to centre spacing 100 um x 100 um, Noise 
Equivalent Power 2.2mW/cm? @ 24Hz. 

- Pyroelectric point detector Spectrum SPH 62 THz - Gentec QS2 THz BL (Ef- 
fective aperture 2 mm x 2 mm, Noise Equivalent Power 4.0 : 107!° W/VHz, Voltage 
Responsivity 140KV/W @ 630nm 5Hz). 

- Pyroelectric point detector Spectrum SPH 31 THz - Gentec QS5 THz MT (Ef- 
fective aperture 5mm Ø, Noise Equivalent Power 1.0: 107°? W/VHz, Voltage Re- 
sponsivity 70KV/W @ 630nm 5Hz). 

- Nanowires detectors (see nanowires paragraph for detailed description). 

THz radiation is not transmitted through glass and, consequently, special optics 
are required; among the various materials currently employed to build THz optics, 
we had at our disposal some Polytetrafluoroethylene (PTFE or Teflon) lenses (quite 
inexpensive and simple to work material with, however, relatively low transmittance 
coefficient) and various Polymethylpentene (TPX) lenses and windows (more ex- 
pensive material with higher transmittance coefficient both in the THz and in the 


66 


Massimiliano Locatelli 


visible range); however, when possible, we decided to use off axis parabolic (OAP) 
mirrors in order to reduce at maximum the beam attenuation. 


2 Experimental setups and results 


Considered the low sensitivity of the array detector at our disposal (Pyrocam III) 
and the limited output power of our QCL, it was not possible to detect collimated 
beams of lateral size comparable to the sensor dimensions and, consequently, use of 
the camera was restricted to focused beam detection and analysis. These detection 
restrictions prevented us from realizing a real time imaging system and forced us to 
implement a single detector system in which the whole sample surface could be in- 
vestigated by means of a motorized scanning system. This kind of imaging configu- 
ration has been widely employed so far in the THz range due to shortage of com- 
mercial detector arrays and, despite its low image acquisition rate, represented a 
convenient solution to test the newly conceived kind of sensors, based on field effect 
transistor nanowires, provided to our laboratory by the Nanoscience Institute of Pi- 
sa. 

The realization of a DH setup without a sufficiently sensitive array detector was 
an even more challenging task. Single detector recording of an interference pattern, 
even if possible in principle [54], would be very critical because of fringe instabilities 
during the scanning process and would require smaller detector elements than ours, 
in order to sample correctly the narrow fringes required in off axis DH. The pur- 
chase of a higher sensitivity array detector was therefore an inevitable decision and a 
microbolometric focal plane array revealed to be the best choice; considered that 
high performance detectors of this kind are not yet widespread, the order required, 
besides great attention and economic effort, long times and it is currently still under 
processing. 


2.1 THz imaging with single point detector 


In figure 4.1 the scheme of the THz imaging setup realized in our laboratory is 
illustrated. We opted for a transmission configuration which could point out THz 
radiation penetration capability. 

The QCL was indium soldered to a copper carrier mounted in the cryorefrig- 
erator where vacuum and low temperatures were obtained. Due to the high angu- 
lar divergence of the emitted beam, the device was placed as close as possible to the 
2 mm TPX cryorefrigerator window. In this configuration the laser diverging out- 
put was collimated by means of a first f/1 OAP mirror with 50.8 mm reflective fo- 
cal length (RFL). The collimated beam was then focused on the sample by means 
of a second f/1 OAP mirror with the same RFL. The sample was fixed, by means of 
a ‘home made’ sample holder frame to two motorized translational stages (LTS300 
by Thorlabs) which could move the sample along the laser beam focus plane. The 
diverging beam coming out from the sample plane was collimated by means of a 
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third f/1 OAP mirror with the same RFL of the previous ones. The beam was final- 
ly focused on a point detector with a fourth, and last, f/1 OAP mirror with 38.1 
mm REL. It should be remarked that this configuration made no use of lenses and 
the beam path had quite contained dimensions in order to reduce at maximum the 
atmospheric attenuation. The motorized stages were remotely controlled by means 
of a ‘home made’ Labview routine so as to obtain a raster scanning of the sample in 
front of the laser beam waist. The attenuated transmitted beam intensity was col- 
lected by the point detector so that image data were acquired one pixel at a time 
and the assembled image of the sample could be recovered at a later stage. Both 
pyroelectric detectors at our disposal and the Pyrocam were used to achieve an op- 
timal alignment. A lock in amplifier (EG&G Instruments 7265 DSP Lock in Am- 
plifier) was used to acquire the point detector response with an improved signal to 
noise ratio. 


parabolic mirror 


parabolic mirror f 50.8mm 


f 50.8mm 


parabolic mirror 
f 38.1mm z] parabolic mirror 
f 50.8mm 


Figure 4.1: Terahertz imaging setup. 


Acquiring with a pyroelectric sensor 

In a first series of measurement, we used, as a point detector, the pyroelectric 
sensor Gentec QS2 THz BL. Working with pyroelectric detectors requires to low- 
frequency modulate the laser beam, e.g. by using a chopper or by modulating la- 
ser driving current. We opted for the second choice by combining the laser driv- 
ing pulses (20 us, at a repetition rate of 10 kHz) with a 50% duty cycle square 
wave at 100 Hz. The detector signal was acquired by means of the lock in ampli- 
fier with a time constant of 20 ms and stored in a *.txt file by means of the same 
Labview routine controlling the motorized stages. One of the two motorized 
stages scanned back and forth the sample surface in front of the focused beam 
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with a fixed speed along the horizontal direction while the lock-in acquired the 
detector signal at a given rate; at each line end the other stage shifted the sample 
in the vertical direction of a fixed quantity. The horizontal and vertical pixel 
pitches of the final image were thus determined, respectively, by the fixed value 
of the vertical step and by the ratio between the horizontal scanning speed and 
the acquisition rate. 

In order to have an estimation of the beam waist symmetry and dimensions 
in the object focus plane, we used, as a first sample, a 250 um pinhole and moved 
it in front of the beam focus with 20 um resolution in both directions covering a 
1 mm x 1 mm surface (graph 4.2). Since the pinhole diameter and the laser beam 
waist had comparable dimensions this method could only give an overestimation 
of the real beam characteristics. The two beam profiles resulted, however, to be 
quite symmetrical and the beam diameter was in good agreement with the ex- 
pected value calculated on the basis of the laser divergence and of the OAP mir- 
ror RFL. 
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Graph 4.2: Beam waist image. 


In order to test the resolution of our imaging system we acquired the image of an 
USAF Target (R74). In this first image (figure 4.2), acquired with a pixel pitch of 
about 100 um x 100 um, we managed to solve lines belonging to group 1 element 4, 
meaning a resolution of 2.828 cycles/mm (about 177 um) comparable to the diffrac- 
tion limited resolution. The unexpected black spots inside some of the target slits are 
acquisition artefacts due to unwanted lock in saturation. 
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Figure 4.2: USAF target THz image. 


We then decided to test our imaging system with a typical sample used in Te- 
rahertz imaging experiments [55], that is a partially dry leaf. The image, acquired 
with a pixel pitch of about 100 um x 100 um, shown in figure 4.3, points out the abil- 
ity of terahertz radiation to pass through opaque organic materials and to highlight 
their water content distribution; imaging with THZ radiation can therefore empha- 
size the vascular bundle system and the internal structure of the leaf. 
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Figure 4.3: Dry leaf THz image. 


This kind of setup can be easily employed in cultural heritage applications as is 
the case, for example, of certain pigments used in painting that are transparent at 
this wavelength, thus allowing the reconstruction of the underlying structures. It is 
thus possible to highlight fractures or detachments in the supporting layers of the 
pictorial work when paints on wood or frescos on masonry surfaces are investigated. 
With our configuration, in particular, we proved the capacity of 2.8 THz radiation to 
penetrate the white lead pigment, (PbCO3).-Pb(OH)2, commonly called biacca; this 
pigment, besides being employed as white painting colour, was sometimes used, 
even by important artists like Michelangelo Buonarroti, to cover pre-existing unde- 
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sired drawings where new drawings were to be realized. We also tested the capacity 
of our system to detect iron gall ink or graphite letters written in old paper (figure 
4.4). With a simple imaging system in the THz range, like the one developed in our 
laboratory, it would thus be possible to bring to light graphite or ink drawing or 
writings intentionally hidden by the artist under a white lead layer. 


Figure 4.4: (Left) Iron gall ink “A” letter. (Right) Graphite “A” letter. 


Acquiring with a nanowire sensor 

In a second session of measurements we had the opportunity to test a new kind 
of THz detectors as point detectors for our imaging system: nanowire detectors. 

Many THz detection system approaches have been developed so far, but a min- 
iaturized room-temperature detector technology, which could be also easily integrat- 
ed in an array configuration, is still largely missing. Semiconducting bolometers and 
hot-electron bolometers can be indeed conveniently coupled with QCLs leading to 
high signal to noise ratio but at the cost of a deep cryogenic cooling and, usually, a 
low dynamic range. On the other hand, room temperature THz detectors, such as 
pyroelectric detectors, Golay cells, and Schottky diodes, are either not very sensitive, 
or extremely slow, or operate well only at IR frequencies lower than 1THz, therefore 
being not ideal for QCL sources. Commercial THz focal plane arrays, currently 
based on microbolometric elements, can be efficiently coupled with QCLs and have 
allowed demonstrating real-time imaging [49;50], but, presently, provide only mod- 
erate sensitivities and response speed. Schottky diodes, the prototypical electronic 
components for room temperature detection of high frequency radiation, exploit the 
nonlinearity and asymmetry of their current-voltage characteristic, to generate, 
through rectification, a continuous signal out of the oscillating incident electromag- 
netic field; this simple operating principle, however, has several drawbacks, related 
to the RC time constant, limiting the bandwidth around 1THz and to the limited 
sensitivity. Considered the large interest in THz imaging applications, the develop- 
ment of a new solid-state technology for fast, room-temperature, terahertz detectors 
integrated in high-speed multi pixel arrays is highly desired. A possible route toward 
such a kind of detectors could be represented by nanostructure devices which, today, 
are at the forefront of studies on future electronics; the growth of semiconductor 
nanowires (NWs), in particular, has recently opened new paths to silicon integration 
of device families such as light-emitting diodes, high-efficiency photovoltaic devices, 
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or high-responsivity photodetectors. NWs can be easily removed from the host sub- 
strate and placed on top of a new functional one for individual contacting, even in 
relatively large numbers, with a simple planar technology suitable for low- 
capacitance circuits. Therefore, in principle, they represent an ideal building block 
for implementing wave detectors that could be well operated into the terahertz, 
thanks to their typical attofarad-order capacitance. Recently, such a novel approach 
for the development of THz detectors by appropriately integrating conventional 
THZ detection geometries with nanowires structures has been proposed [56;57]. 
Specifically, semiconductor InAs nanowires grown bottom-up by catalytic process 
through vapour-phase epitaxy has been used as active element for the development 
of THZ detectors based on a FET configuration. The detection mechanism exploits 
the nonlinearity of the transfer characteristics: the terahertz radiation field is fed at 
the gate-source electrodes through wide band antennas, and the rectified signal is 
then read at the output in the form of a DC drain voltage. The terahertz detection 
principle in a FET was first explained by the Dyakonov-Shur plasma wave theory 
[58]. When terahertz radiation is coupled to the FET, between gate and source ter- 
minals, the ac terahertz voltage modulates simultaneously the carrier density and the 
carrier drift velocity. As a result, the terahertz signal is rectified and leads to a DC 
signal AV between source and drain terminals proportional to the received power. 
The value of this voltage (or current, depending on the read-out circuit) depends on 
the carrier density in the channel, which may be controlled by the gate voltage. The 
device operates as a square law detector in which the largest non-linearity, leading to 
the highest responsivity, is achieved around the channel pinch-off. In addition, some 
asymmetry between source and drain is needed to induce AV; this can originate 
from the difference in the source and drain boundary conditions due to some para- 
sitic capacitance, but usually stems from the asymmetry in feeding the incoming ra- 
diation, which can be achieved either by using a special antenna or by an asymmetric 
design of source and drain contact pads. Finally, the asymmetry can naturally arise if 
a dc current is passed between source and drain, creating a depletion of the electron 
density on the drain side of the channel. In particular, to increase the asymmetry, 
and therefore the responsivity, of the nanowire detectors provided to our laboratory, 
low shunt-capacitance antennas to funnel the radiation into the strongly sub wave- 
length detecting elements, were designed as contact pads. Antenna coupling ensures 
selective responsivity to both the spatial mode and the polarization of the incoming 
radiation. These very sensitive sensors, as every FET device, are quite delicate and 
need particularly careful handling especially from the point of view of accidental 
electrostatic discharges. 

Following the standard characterization procedures performed in the papers de- 
scribing these devices, we decided to test the electrical and responsivity properties of 
these nanowires. By means of X,Y,Z micro movements stages, we positioned the 
nanowire detector in the first focus of the beam path, where higher power was at- 
tained. We drove the laser with a 5.04 kHz and 30% duty cycle square wave current 
with amplitude of 0.75 A at a temperature of 25 K (when working with NWs burst- 
ing is no longer required). By means of the stages we positioned the NW so as to op- 
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timize the radiation-antenna matching and obtained a maximum signal of 350 uV at 
the drain electrode while applying a 7 V DC voltage to the gate terminal. The signal 
was then amplified 10X and sent to a real time signal analyser (Tektronix RSA 
51106A) so that we could extract the power spectral density in a 200 Hz window and 
with 1 Hz of resolution bandwidth in three different conditions: laser on, laser off, 
laser on but blocked beam. 
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Graph 4.3: Power density spectrum of the signal with laser on (red), with laser off (blue) and 
with laser on but blocked beam (black). 


As shown in graph 4.3 we had about 77 db/Hz of SNR even if a portion of the 
signal at the working frequency was probably coming from the laser electronic cir- 
cuits. 

We then acquired the frequency response of the detector in a broader frequency 
window (graph 4.4) gradually changing the frequency of the laser driving current. 


77dB/Hz 
—— QCL signal 
Detector noise 
—— 20k resistor noise 


Power spectral density (dBm/Hz) 


10 100 1000 10000 100000 1000000 


Frequency (Hz) 


Graph 4.4: (Blue) Nanowire frequency response. (Red) Power spectral density of the detector 
noise. (Black) Power spectral density of a 20 KQ resistor. 
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In the same graph it is also represented the power spectral density of the detector 
noise and the power spectral density of the Johnson noise at the terminals of a 20 kQ 
resistor in the same frequency bandwidth. Coherently with the results presented in 
[57], the frequency response of the detector remained constant within 3 db up to 
about 5 kHz. In this frequency region the detector noise spectrum seemed to be cen- 
tered around about 115 dbm/Hz as if the detector resistance were of the order of 120 
kQ. 


We then focused our attention to the detector response behaviour for different 
gate voltage values. To investigate this relation we used a waveform generator to ap- 
ply a triangular wave at 1.4 Hz going from -8 V up to 8 V to the gate terminal. As 
shown in graph 4.5, the maximum nanowire response was obtained for a positive 
value of the gate potential V, (about 6 V) but the curve showed clear hysteresis be- 
haviour as confirmed in [56;57]. 
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Graph 4.5:.Nanowire signal vs gate potential. 


In the last characterization experiment we measured the source-drain I-V char- 
acteristics and recovered the resistance of the source-drain channel for different val- 
ues of Vy. To perform this measurement we applied a constant voltage to the drain 
terminal in series with a 5 MQ resistor, and measured the current intensity at the 
output of the source terminal by means of a current amplifier converting the current 
into a voltage signal with an amplification factor of 10° V/A. The latter signal was 
sent to the oscilloscope together with the triangular wave at 10 mHz going from -8 V 
up to 8 V applied to the gate terminal. The results are illustrated in graphs 4.6 and 
4.7. 
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Graph 4.6: Source-drain I-V characteristic. 
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Graph 4.7: Source-drain channel resistance Ras vs gate voltage Vs. 


Finally we put back the detector on the second beam focus and started acquiring 
images with our nanowire THz imaging setup. We first acquired an image of the 
usual USAF target to verify the image resolution. As shown in figure 4.5 the result 
was comparable with the result obtained using the pyroelectric sensor (here, actually, 
we had a little bit higher resolution because of a better beam focusing and arrived to 
solve completely the elements belonging to group 1 element 6, corresponding to a 


resolution of about 140 um). 
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Figure 4.5: USAF target THz image using nanowire detector. 


Good results were obtained using a dry leaf as a sample, acquired with the usual 
pixel pitch of about 100 um x 100 um (figure 4.6). 
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Figure 4.6: Dry leaf THz image using nanowire detector. 


The striped background clearly observable in both images is due to the instabil- 
ity of nanowire response on long times (of the order of 5% in10 s), an effect that is 
still under investigation. 


2.2 Toward Terahertz digital holography 


Terahertz Digital Holography is still an almost unexplored research field. Few 
results, in fact, have been obtained in this direction. The first experiment on Far IR 
DH dates back to 2006 when an off-axis transmission digital holography setup at 
millimetre wavelenghts was realized using a 100 GHz Gun diode oscillator as a 
source and a spatially scanned Schottky diode as detector [59]. In 2008 computer 
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simulations and an experimental realization of THz transmission digital holograms 
were realized with a THz TD technique [60]. Further transmission holography con- 
figurations were realized in 2011 using a 50 mW Far Infrared gas laser at 2.52 THZ 
and a pyroelectric camera [61] and using a frequency tunable continuous wave THz 
source emitting at 0.7 THz and a spatially scanned schottky diode detector [54]. To 
date, however, no speckle holographic configuration has ever been realized and QCL 
THZ sources have never been used to this purpose. Our next step is to develop such 
a configuration using our QCL and a microbolometric camera: considered the low 
output power of our laser, a high sensitivity array is needed and therefore we decided 
to buy a focal plane array microbolometric camera, IRXCAM INO 160 THz [62;63]; 
this camera is composed of 160 x 120 VaO square pixels with 52 um pixel pitch and 
is optimized for the THz region detection in order to reach a NEP of 100 pW. 
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Figure 4.7: Simple THz holographic setup. 


Our starting point will be the investigation of small scattering objects with a very 
simple and power saving configuration like the one illustrated in figure 4.7. In this 
configuration the laser beam is collimated by an OAP mirror and then it is inter- 
cepted both by the small object under investigation and by an adjustable mirror; the 
portion of the beam impinging on the sample, constitutes, as usual, the object beam 
while the portion of the beam reflected by the mirror constitutes the reference beam. 
With this configuration it is possible to obtain the desired fringe spacing simply 
changing the distance between the microbolometric camera and the sample and ad- 
justing the reference mirror inclination. This simple configuration will be used to 
demonstrate the feasibility of the techniques but in a second step we will try to trans- 
late all the configurations developed in the field of Mid IR DH to the THz domain in 
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order to take full advantage of the stability and of the increased field of view offered 
by a THz radiation DH setup. 
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Chapter 5 
Conclusions 


During this research project development I had to face various arguments in the 
laser optics and imaging field. I had the opportunity to work on well-known lasers 
sources (CO; lasers) and detectors (pyroelectric cameras and detectors, microbolo- 
metric cameras) as well as on frontier sources (Quantum Cascade lasers) and inno- 
vative detectors (nanowires). I learnt to exploit the different characteristics of these 
specific coherent radiation sources and sensors to implement new and promising 
imaging techniques. 

Nowadays we are witnessing a continuously growing interest in the Mid IR re- 
gion and, even more, in the Far IR region, especially for imaging applications. Mid 
IR imaging systems are widely employed in the military and security field but also 
for energy saving and medical investigations and are, consequently, in continuous 
development. Even greater expectation is placed on imaging techniques in the Te- 
rahertz region, because of its well-known capacity to pass undisturbed through many 
common materials and to provide spectroscopic information about various strategic 
materials like explosives or drugs. In this scenario DH, a quite recent and peculiar 
imaging technique, is proving to be mature enough to play a key role among the 
numerous imaging techniques both in the thermal IR and in the Terahertz IR, espe- 
cially for some specific applications. As we have seen, IRDH at 10.6 um has paved 
the way to large size object investigation, especially for non-destructive testing and, 
thanks to its ability to see through smoke, dust particles and flames, it could be em- 
ployed, in the future, in military and safety applications. Mid IRDH could play a very 
important role also in the research activity about 3D display: holographic television, 
commonly regarded as the holy grail of holography, is one of the most promising 
and challenging developments for the future display market; only holography can 
provide all depth information necessary to obtain the reconstruction of natural- 
looking 3D scenes; IRDH, thanks to its capacity to record whole wavefront infor- 
mation of large size samples in daylight and not isolated conditions, and thanks to 
the possible developments introduced by the progresses of spatial light modulators, 
may have an important role in this challenge. A desirable step to increase Mid IRDH 
could be accomplished if more compact sources, like QCLs, should replace the in- 
credibly versatile but equally bulky CO: lasers. An important purpose for the future 
is to pursue these results and to repeat all the successful experiments carried out in 
the Mid IR range in a miniaturized, and therefore more appealing and usable, ver- 
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sion. All the results obtained in the Mid IR DH activity could be fruitfully extended 
to the THz domain with improved results: an even lower vibration sensitivity and 
larger field of view should be obtained using a 10 times longer wavelenght than the 
one we have used up to now; further, the repeatedly highlighted capability of THz 
radiation to penetrate many materials impenetrable to visible light, combined with 
the capacity of DH to reconstruct the 3D images, would represent an invaluable tool 
to analyse internal structure with 3D perspective, otherwise not achievable. All these 
applications will receive an important impulse from the development of new and 
more sensitive and fast IR detectors like, hopefully, the nanowires FET tested with 
our imaging THz setup. 
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